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1.  INTRODUCTION 


This  is  the  final  report  on  Air  Force  Phillips  Laboratory/Geophysics  Directorate 
Contract  F19628-92-C-0016.  Work  emphasized  the  conversion  of  a  suite  of  CPI  research 
codes  into  a  new  suite  using  modem  programming  standards  involving  redesign,  software 
implementation,  testing,  and  validation.  The  original  statement-of-work  (SOW)  called 
for  a  new  suite  of  codes  that  would  model  spectral  radiances  for  dayglow,  nightglow,  and 
electron  aurora.  Due  to  inteiruptions  in  funding  and  associated  uncertainty  about  the 
needed  support  to  complete  the  original  SOW,  the  program  was  de-scoped  to  exclude 
auroral  modules.  The  final  products  are  a  suite  of  codes  that  model  dayglow  and 
nightglow  along  with  documentation.  The  latter  includes  a  User’s  manual  and  two 
science  papers.  The  papers  follow  this  introduction.  We  have  just  received  word  that  the 
second  paper  titled  “New  Survey  of  Electron  Impact  Cross  Sections  for  Photoelectron  and 
Auroral  Electron  Energy  Loss  Calculations”  as  been  accepted  for  publication  in  J.  Phys. 
Chem.,  Ref.  Data. 

The  delivered  software  is  to  be  part  of  a  larger  set  of  models  called  AURIC  which 
stands  for  Atmospheric  Ultraviolet  Radiance  Integrated  Code.  The  term  Integrated  refers 
to  original  plans  to  develop  a  user  interface  behind  which  would  lie  the  delivered 
software  along  with  MODTRAN  and  associated  software.  The  combined  capability 
would  allow  users  to  calculate  spectral  radiances  from  the  upper  atmosphere  (airglow), 
from  the  middle  and  lower  atmospheres  (Rayleigh  scattering  of  sunlight  and  moonlight), 
as  well  combinations  of  these  emissions.  Plans  are  imcertain  at  this  time  with  regard  to 
this  integration. 

The  delivered  software  offers  a  stand-alone  capability  to  calculate  airglow  spectral 
radiances  from  100  to  900  nm.  Viewing  is  from  satellite  altitudes  and  includes  the  limb 
and  hard  disk.  The  software  was  redesigned  in  Fortran  and  is  available  for  running  on 
Unix,  VMS,  and  PC  systems.  All  input  and  output  files  are  formatted  with  standard 
formats  wherever  possible.  The  suite  of  codes  is  comprised  of  27  members  containing 
approximately  35,000  lines  of  coding.  It  provides  the  user  with  the  capability  to 
characterize  optical  backgrounds  and  conduct  detailed  scientific  studies.  Studies  include 
interpretation  of  space-based  radiance  measurements  and  theoretical  investigations  of  a 
variety  of  optical  features,  neutral  densities,  and  ion  densities. 

Table  1  lists  atomic  dayglow  and  nightglow  features  presently  contained  in 
AURIC.  Similar  information  for  band  systems  is  given  in  Table  2.  Emission  for  some  of 
the  dayglow  and  nightglow  features  arises  from  metastable  states.  Densities  of  these 
states  are  calculated  by  AURIC’s  photochemistry  and  nightglow  chemistry  models. 
Table  3  lists  the  dayglow  metastables  which  are  part  of  a  larger  list  that  includes  ground 
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state  ion  species  and  N(  S). 


The  original  fimding  for  this  project  came  from  the  Ballistic  Missile  Defense 
Office  with  the  intent  of  providing  analysts  with  a  tool  for  accurately  calculating  optical 
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backgrounds  as  seen  from  space.  The  delivered  software  achieves  this  goal  for  two  of  the 
three  upper  atmospheric  emission  components  (dayglow  and  nightglow).  Since  first- 
principles  rather  than  empirical  approaches  have  been  used  to  characterized  the  radiances, 
the  software  has  a  broader  use.  Quantitative  analysis  of  satellite  and  rocket  data  may  be 
undertaken  for  optical  emissions  as  well  as  ion  and  neutral  densities.  The  software  may 
also  be  used  in  the  development  of  remote  sensing  algorithms.  An  example  is  an 
algorithm  developed  by  Strickland  et  al.  (J.  Geophys.  Res.,  100,  12,217,  1995)  for 
deriving  the  abundance  of  atomic  oxygen  and  an  integrated  value  of  the  solar  EUV  flux 
from  N2  Lyman-Birge-Hopfield  and  01 135.6  nm  disk  viewing  dayglow  data. 


Table  1 .  Atomic  lines,  transitions,  type  of  airglow  (day/night)  modeled,  and  identification  of  those  features 
affected  by  scattering  or  requiring  chemistry  modeling  (metastables) 


Species 

Transition 

Wavelength 

Region 

Comment 

(nm) 

H 

y->'s 

121.6 

Global 

Requires  multi-dimensional  photon 

transport  model 

N 

149.3 

Day 

N 

■?  ft 

174.3 

ti 

N 

p^p 

346.6 

M 

Metastable 

N 

P  s 

520.0 

II 

Metastable 

D  S 

N+ 

D  P 

108.5 

Day 

N+ 

214.3 

11 

Metastable 

0 

S  P 

130.4 

Day/Night 

Requires  photon  transport  model 

0 

135.6 

II 

Requires  photon  transport  model 

0 

b  ^  r 

297.2 

II 

Metastable 

0 

’s-^'p 

557.7 

M 

Metastable 

0 

's-^'d 

630.0 

II 

Metastable 

0 

0 

'd-^'p 

777.4 

844.6 

11 

ft 

p-^  s 

P“^  s 

0+ 

p-^  s 

247.0 

Day 

Metastable 

0+ 

2_o _ v4_o 

D  S 

372.7 

II 

Metastable 

732.0 

Metastable 

0+ 

3„o_,2^o 

P  D 

-2- 


A  conscientious  effort  was  made  to  test  and  validate  the  software.  Testing 
included  extensive  code  reviews  and  cross  platform  verification  of  outputs.  Validation 
included  comparisons  with  results  from  the  original  research  codes  (which  themselves 
have  been  validated  over  the  years  through  model/data  comparisons)  and  comparisons  to 
selected  emission  and  composition  data. 

Considerable  effort  was  directed  to  the  User’s  manual.  Every  I/O  file  and  code 
including  a  user  interface  is  discussed.  Detailed  information  is  provided  on  how  to  run 
the  codes  on  each  of  the  three  types  of  platforms.  Detailed  I/O  information  is  also  given 
in  the  form  of  plots  of  key  input  parameters  and  calculated  quantities  such  as  the 
photoelectron  source  function,  the  photoelectron  flux,  volume  production  rates,  limb 
intensity  profiles,  spectral  radiances,  and  densities  of  ions  and  neutrals. 


Table  2.  Band  systems,  transitions,  type  of  airglow,  and  wavelength  regions  covered. 


Species 

Transition 

Band  system 

Region 

Wavelength  range  (nm) 

N2 

+ 

a  rio  X  Za 

0  & 

Lyman-Birge-Hopfield 

Day 

125-240 

N2 

a'Xu"  x'Zg" 

Vegard-Kaplan 

t! 

150-690 

N2 

C’Hu  B^Rg 

Second  Positive 

H 

280-460 

N2 

B  ria  A  Zy 
& 

First  Positive 

M 

600  to  beyond  1000 

N2+ 

2  +  2  + 

B  Zy  — >  X  Zg 

First  Negative 

Day 

360  -  520 

N2'" 

a'Ru  ^  x'Zg* 

Meinel 

fi 

500  to  beyond  1000 

NO 

aV  x'n 

y 

Day 

200-310 

NO 

x"n 

5 

tt 

190-250 

NO 

dV  x’n 

8 

180-260 

02 

b'Zo"  x'z„- 
&  & 

Atmospheric 

Night 

680  -  1000 

O2 

/j  i 

A  Ay  — ^  a  Ag 

Chamberlain 

»» 

300  -  440 

O2 

A  Z  — ^  X  Zcr 

Herzberg  I 

250-480 

O2 

c'Zy--^  x'Zg- 

Herzberg  II 

It 

260  -  470 

As  a  concluding  introductory  comment,  much  more  information  on  the  overall 
capability  of  the  delivered  software  may  be  seen  in  the  first  of  the  following  two  papers. 


We  would  like  to  acknowledge  Dr.  Robert  Huffman  of  AFPL/GP  who  has  been  the 
COTR  of  this  project  during  its  entire  lifetime.  His  support  and  comments  have 
strengthen  the  program  and  helped  lead  to  completion  of  the  revised  SOW. 


Table  3.  Species  treated  by  the  photochemistry  module  and  identification  of  dominant  production 

mechanisms 


Species 

Dominant  Production  Mechanisms 

Slow  reacting  species  (time  integrated) 

NO* 

Charge  transfer 

O3 

photoionization,  charge  transfer 

0* 

e  impact,  photoionization 

N('S) 

N(“D)  quenching  by  0 

Fast  reacting  species 

N* 

e  impact,  photoionization 

K 

e  impact,  photoionization 

0*CD) 

e  impact,  photoionization 

+  2 

0(P) 

e  impact,  photoionization 

N('d) 

e  impact,  photodissociation 

N('P) 

e  impact,  photodissociation 

o('d) 

e  impact,  O^^recombination 

0('s) 

e  impact,  quenching  of  N^CA)  state 

N,(A) 

e  impact 
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APPENDIX  A:  AURIC  THERMOSPHERIC-IONOSPHERIC  CHEMISTRY 

RATE  COEFFICIENTS 
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Modeling  of  thermospheric-ionospheric  chemistry  within  AURIC  is  geared  primarily  towards  evaluation 
of  volume  excitation  rates  of  prominent  metcistable  species  and  associated  /olume  emission  rates  from  radiative 
relaxation.  The  chemical  reaction  scheme  is  fairly  complete  with  regard  to  N,  0  and  0"^  metastables;  of  the 
molecular  species,  the  only  metastable  explicitly  followed  is  N2(A).  The  reaction  rate  coefficients  and  pertinent 
references  are  given  in  Tables  A1-A4.  Alternate  recent  discussions  of  ionospheric-thermospheric  chemistry 
within  the  context  of  aeronomic  data  analyses  are  provided  by  Link  et  al  [1983],  Fox  and  Dalgarno  [1985], 
Link  and  Cogger  [1988],  Solomon  ei  ai  [1988],  Buonsanto  ei  al  [1992],  and  Torr  et  al  [1995]. 

By-and-large,  reaction  rate  coefficient  expressions  have  been  taken  from  laboratory  measurements  where 
available,  even  in  those  cases  where  past  aeronomic  data  analyses  may  have  suggested  reassessment  of  these 
rates.  Pertinent  rate  coefficients  for  some  reactions  (i.e.,  rate  coefficients  for  genuine  thermal  distributions  at 
temperatures  250  K  ^  T  ^  1500  K),  however,  have  not  been  experimentally  determined;  in  these  cases,  rates 
derived  from  aeronomic  data  analyses  have  been  used  where  available  {e.g.,  Consiantinides  ei  al  [1979]).  In 
several  cases,  reaction  rate  coefficients  that  have  been  treated  as  “standard”  in  aeronomic  modeling  have  been 
questioned  in  recent  laboratory  investigations.  Examples  are:  Flesch  and  Ng  [1991]  vs  Johnsen  and  Biondi 
[1980a,b],  Schultz  and  Armenirout  [1991]  vs  McFarland  et  al  [1973],  and  Burley  et  al  [1987]  vs  Chen  et  al 
[1978].  We  have  retained  the  “standard”  rate  coefficients  in  these  cases,  in  part  because  the  more  recent 
laboratory  results  are  not  convenient  and  in  part  because  the  “standard”  rates  have  been  adequate  for  previous 
aeronomic  applications. 

Since  the  vibrational  distributions  of  several  molecular  species  are  not  well  characterized  Nj,  Oj, 
NO“*“,  N2(A)),  the  rates  adopted  here  are  biased  towards  the  relevant  ground  vibrational  level.  For  example, 
the  dissociative  recombination  rate  and  branching  ratios  for  production  of  product  metastable  species  for 

0+  +  e-  0{^D)^0{^P) 

0(^D)  +  0(^T>) 

0(^1))  + 0(^5) 

are  from  Meyr  and  Biondi  [1969]  with  the  branching  ratios  for  the  product  channels  taken  from  Guherman 
[1988]  for  the  Oj  =  0  state.  There  is  no  tracking  of  translationally  excited  species. 

Miscellaneous  Thermospheric-Ionospheric  Chemistry  References 

Buonsanto,  M.  J.,  S.  C.  Solomon,  and  W.  K.  Tobiska,  Comparison  of  measured  and  modeled  solar  EUV  flux 
and  its  effect  on  the  E-Fl  region  ionosphere,  J.  Geophys.  Res.,  97,  10513-10524,  1992. 

Fox,  J.  L.,  and  A.  Dalgarno,  The  vibrational  distribution  of  Nj  in  the  terrestrial  ionosphere,  J.  Geophys.  Res., 
90,  7557-7567,  1985. 

Link,  R.,  and  L.  L.  Cogger,  A  reexamination  of  the  Ol  6300-A  nightglow,  J.  Geophys.  Res.,  93,  9883-9892, 
1988. 

Solomon,  S.  C.,  P,  B.  Hays,  and  V.  J.  Abreu,  The  auroral  6300  A  emission:  Observations  and  modeling,  J. 
Geophys.  Res.,  93,  9867-9882,  1988. 

Torr,  M.  R.,  D.  G.  Torr,  T.  Chang,  P.  Richards,  W.  Swift,  and  N.  Li,  Thermospheric  nitric  oxide  from  the 
ATLAS  1  and  Spacelab  1  missions,  J.  Geophys.  Res.,  100,  17389-17413,  1995. 
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Table  Al:  AURIC  Chemical  Reactions 


Reaction 


Rate  Coefficient  (cm^  s"^) 


References 


N++02 

^  Oj  +  N2 

fci  =  5.0  X  10"“  7--°  ® 

McFarland  et  al.  [1973] 
Lindinger  ei  al  [1974] 
Smith  ei  al  [1978] 

N+  +  NO 

N0+  +  N2 

ki  =  3.3  X  10-1“ 

Fehsenfeld  ei  al  [1970] 

N+  +  0 

— >  products 

ks  =  1.4  X  10-1“  Ti  <  1500“K 

=  5.2  X  10-11  Ti“-2,Ti>  1500”K 

McFarland  ei  al  [1974] 
Knuisen  ei  al  [1988] 

N0+  +N(2U'’) 

kaa  =(l-0.07T-“-2i)-fc3 

0+  +  N2 

k3b  =0.07Ti“-2i.l;3 

Oj  +  NO 

-  N0+  +  O2 

ki  =  4.4  X  10-1“ 

Lindinger  ei  al  [1975] 
Glosik  ei  al  [1978] 
Lindinger  ei  al  [1974] 

0J  +  N(45‘>) 

—  N0+  +  0 

ks  =  1.2  X  10-1“ 

Fehsenfeld  [1977] 

O+  +  N2 

N0+  +  N(45‘>) 

ks  =  2.78  X  10-1“  exp(2.07/7',  —  0.608/t-“) 

Chen  ei  al  [1978] 

Smith  ei  al  [1978] 
Lindinger  ei  al  [1974] 
McFarland  ei  al  [1973] 

O+  +  O2 

-  0+  +  0 

jfcr  =  3.25  X  10-11  -  1.30  x  lO-n  r. 

+  1.60  X  10-1“  7-“ 

Chen  ei  al  [1978] 

Smith  ei  al  [1978] 
Rakshii  ei  al  [1978] 
Lindinger  ei  al  [1974] 
McFarland  ei  al  [1973] 

0+  +  N0 

N0+  +  0 

ks  =  8.0  X  10-1“ 

Craham  et  al  [1975] 
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N+  +  02 

-H.  products 

kg  =6.1x10-1° 

Adams  ei  al  [1980] 

Rakshii  ei  al  [1978] 
Lindinger  ei  al  [1974] 
McFarland  ei  al  [1973] 

N0+  +  0 

N0++0(1D) 

kga  =  0.3  •  0.43  •  ikg 

kgb  =  0.7  ■  0.43  ■  kg 

Adams  ei  al  [1980] 

Langford  ei  al  [1986] 

0j  +  N('>5‘>) 

Oj  +  N(2d‘>) 

kgc  =  0.7  •  0.51  •  kg 

kgd  —  0.3  •  0.51  •  kg 

Adams  ei  al  [1980] 

O’Keefe  ei  al  [1986] 

0+  +  NO 

kge  =  0.06  •  kg 

Adams  ei  al  [1980] 

N+  +  N0 

products 

N0+  +  N{^5®) 

-  N+  +  0 

^tio  =  5.3  X  10“i° 

kxQd  —  0.85  ■  1*10 

^106  “  0.15  ■  kiQ 

Adams  ei  al,  [1980] 

Fahey  ei  al  [1981] 

N++0 

N(45®)  +  0+ 

ifcii  =  2.2  X  10-12 

Consianiinides  ei  al  [1979] 
Baies  [1989] 

0+(2po)^N2 

Nj  +  0 

ki2  =  4.8  X  10-1° 

Rusch  ei  al  [1977] 

0+(2p'>)  +  02 

0++0 

*13  =  4.8  X  10-1° 

Link  ei  al  [1983] 

0+(2po)  +  0 

0+(2d®)  +  0 

*14  =  5.2  X  10-11 

Rusch  ei  al  [1977] 

Rusch  and  Gerard  [1980] 

0+(2£)°)+N2 

-  N++0 

*15  =8.0x10-1° 

Johnsen  and  Biondi  [1980a, b] 
Rowe  ei  al  [1980] 

0+(2D°)+02 

0+  +  0 

*16  =  7.0  X  10-1° 

Johnsen  and  Biondi  [1980a, b] 
Rowe  ei  al  [1980] 

0+(2£)o)  +  NO 

^  N0+  +  0 

*17  =  1.2  X  10-° 

Glosik  ei  al  [1978] 

0+(2D')  +  0 

0+  +  0 

*18  =  5.0  X  10-12 

Abdou  ei  al  [1984] 
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N(2P“)  +  N2 

N(45°)  +  N2 

kiQ  =  2.0  X  10-1® 

Schofield  [1979] 

N(2p»)+02 

N0  +  0(i£>) 

ho  =  2.2  X  10-12 

Piper  [1993a] 

Schofield  [1979] 

N(2p»)  +  N0 

N(2£>'')  +  N0 

^21  =  3.0  X  10-11 

Schofield  [1979] 

N(2p'>)  +  0 

N(2d<>)  +  0 

k22  =  1.7  X  10-11 

Piper  [1993a] 

Bates  [1989] 

N(2po)^N(^5‘’) 

N(2d*)  +  N(^5'’) 

*23  =  6.0  X  10-1® 

Young  and  Dunn  [1975] 
Taghipour  and  Brennen  [1979] 

N(2D'>)  +  N2 

^  N('‘5'’)  +  N2 

*24  =  1.0  X  10-1®  exp(-510/rn) 

Slanger  and  Black  [1976] 

N(2Z)»)+02 

NO  +  0 

*25a  =  5.9  X  10-12 

Fell  et  al  [1990] 

Schofield  [1979] 

NO  +  0(  ^D) 

*256  =  6.0  X  10-1® 

Link  [1983] 

Link  and  Swaminathan  [1992] 

N(2D'’)  +  N0 

N2  +  O 

*26  =  6.7  X  10-11 

Fell  ei  al  [1990] 

N(2Z)°)  +  0 

N(^5‘>)  +  0 

*27a  =  6.9  X  10-1® 

Fell  ei  al  [1990] 

Piper  [1989b] 

NO++e- 

*276  =  2.5  X  10-1®  (2205  +  Tn)  exp(-4410/T„) 

Ringer  and  Gentry  [1979] 

N(2D'>)  +  0+ 

N+  +  0 

*28  =  1.3  X  10-1° 

Const antinides  et  al  [1979] 
Bates  [1989] 

N(^5‘>)  +  02 

NO  +  O 

*29  =  5.48  X  10-12  exp(-3221/Tn) 

Becker  ei  al  [1969] 

N(45“)  +  N0 

N2  +  O 

*30  =  2.2  X  10-11  exp(160/r„),  Tn  <  400'’K 
=  3.3  X  10-11,  r„  >4oo'>K 

Wennberg  ei  al  [1994] 

Lee  ei  al  [1978] 

N(V)  +  0 

NO  +  Ai/(NO  IR) 

*31  =  3.33  X  10-1®  r-1/2  _  o.567rn"^^^) 

Sun  and  Dalgarno  [1992] 
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0(^5) +  02 


products 


kz2  =  2.32  X  10-12 


Capeianakis  ei  ah  [1993] 


— >  O  +  O2 

— >  0(  +  O2 

0(^5)  + NO 

— »■  products 

^  O  +  NO 

0(  ^D)  +  NO 

0(i5)  +  0 

^  0  +  0 

0(1D)  +  N2 

^  O  +  N2 

0(ll>)  +  02 

— ►  0  +  O2 

0(  ^D)  +  NO 

O  +  NO 

0(  ^D)  +  0 

-.0  +  0 

N2(j4)  +  O2 

— *■  N2  +  O2 

N2(^)  +  N0 

Nz  +  NO 

N2(^)  +  0 

— »■  products 

-  N2  +  0(15) 

N2  +  O 

N2(A)  +  N(^5“) 

— products 

^  N2  +  N(2p<>) 

N2  +  N(2Z?®) 

kz2a  =  0.69  •  kz2 
kz2b  =  0.31  •  kz2 

kss  =  8.0  X  10-11 
^33a  =  0.36  •  ^33 
^336  =  0.64  •  Ar33 

Ar34  ==  2.0  X  10”1^ 

kss  =  2.0  X  10-11  exp(108/Tn) 

ks6.  =  2.9  X  IQ-ii  exp(67.4/Tn) 

*37  =  1.5  X  10-1° 

*38  =  7,0  X  10-12 

*39  =2.5x10-12 

*40  =  5.6  X  10-11 

*41  =2.8x10-11 
^41a  ”  0.75  ’  *41 
*416  =  0.25  •  *41 

*42  =4.0x10-11 
^42a  =  0.9  •  *42 
*426  =  0.1  *  *42 


Slanger  and  Black  [1978b] 

Filseih  ei  al  [1970] 

Slanger  and  Black  [1978a] 

Krauss  and  Neumann  [1975] 

Sireit  ei  al  [1976] 

Sireii  ei  al  [1976] 

Young  ei  al  [1968] 

Ahreu  ei  al  [1986] 

Yee  ei  al  [1990] 

Thomas  and  Kaufman  [1985] 
Thomas  ei  al  [1987] 

Piper  ei  al  [1981a] 

Thomas  ei  al  [1987] 

Piper  ei  al  [1981b] 

Thomas  and  Kaufman  [1985] 
Piper  [1982] 

De  Souza  ei  al  [1985] 

Piper  [1989a] 
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Notes: 

Ti  =  ^5  where  Ti  is  ion  temperature.  Tn  is  neutral  temperature. 

-  ki:  Alge  and  Lindinger  [1981]  indicate  rate  is  insensitive  to  v'  value.  See  Schultz  and  Armentrout  [1991]  for 
recent  laboratory  results  suggesting  a  need  for  revision  in  adopted  rate. 
ksa'  Frederick  and  Rusch  [1977]  suggest  N(^Z)®)  production  must  be  maximum  allowable. 
kA:  Lindinger  and  Ferguson  [1983]  indicate  rate  is  insensitive  to  v*  value. 
ks:  Possible  0(^5)  production  ignored. 

ke:  See  Burley  et  al  [1987]  for  recent  laboratory  results  suggesting  a  need  for  revision  in  adopted  rate. 

^9:  Small  branching  ratio  for  0(  ^5)  production  ignored. 

fcio:  Possible  production  of  O"^  ignored. 

ki2,ki^:  Rusch  et  al.  [1977]  inferred  loss  rates  only. 

See  Flesch  and  Ng  [1991]  for  measured  cross  sections  at  above-thermal  energies  indicating  a  need  to 
reevaluate  the  commonly  accepted  rate  coefficient  values.  Possible  NO"^  [Glosik  et  ah,  1978]  and  N"*” 
production  channels  ignored. 

^13:  Assigned  by  analogy  with  +  N2  [Rusch  et  ai,  1977], 

^17:  Glosik  et  al.  [1978]  did  not  distinguish  between  and  here  taking  their  rate  as  applying 

to 

kis:  Estimated  upper  limit.  Possible  O(^Z^)  production  ignored. 

^19:  Simple  quenching  assumed. 

^20-  Possible  0(  ^5)  production  ignored. 
k2i:  Simple  quenching  assumed. 

k22‘  Following  arguments  of  Bates  [1989]  regarding  allowed  product  channels. 
k25b-  Estimated  upper  limit. 

^26*  Possible  0(^Z>)  and  0(^5)  production  ignored. 
k27a‘  Possible  0{^D)  production  ignored. 
k^A-  Simple  quenching  assumed. 

^36:  Production  of  02(^)  currently  ignored. 

Arss:  Estimate  consistent  with  Baluja  and  Zeippen  [1988]  transition  probabilities  and  Guberman  [1988]  v*'  =  0 
0(  ^D)  branching  ratio. 

^39~^42:  Assuming  N2(A)  in  the  v"  =  0  state;  dayside  mean  vibrational  population  not  known. 

kss-  Simple  quenching  assumed;  Golde  and  Moyle  [1985]  indicate  N2  +  O  -f  O  to  be  a  major  product  channel. 

k4o:  Simple  quenching  assumed. 

kA2'  Branching  ratios  assigned  somewhat  arbitrarily;  Piper  [1989a]  indicates  N(^P°)  production  channel  is 
dominant. 
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Table  A2;  AURIC  Molecular  Ion  Dissociative  Recombination  Rates 


Reaction 

Rate  Coefficient  (cm®  s“^) 

References 

Ot  +  e- 

products 

ai  =  1.95  X  10"'^  r-®  ’’,  <  1200°K 

=  1.61  X  10-^  Te  >  1200'’K 

Mehr  and  Biondi  [1969] 

Alge  et  ai  [1983] 

^  0(ii?)  +  0(3p) 

Qia  =  0.609  •  ai 

Guberman  [1988] 

O0D)  +  O{^D) 

au  =  0.389-01 

Oi^S)  +  OCD) 

aic  =  0.002  •  ai 

N0+  +  e- 

products 

02  =  4.20  X  10“’^ 

Dulaney  et  ai  [1987] 
Davidson  and  Hobson  [1987] 
Alge  et  ai  [1983] 

^  N(^5°)  +  0 

Ci2a  =  0.24  •  02 

Kley  et  ai  [1977] 

N(2£)°)  +  0 

02i  =  0.76  •  02 

N++e- 

products 

03  =  1.80  X  10-"^  r-o  ®® 

Mehr  and  Biondi  [1969] 

N(2£)°)  +  N(^5“) 

030  =  0.88  •  03 

Guberman  [1991] 

N(2£)°)  +  N(2d°) 

03J  =  0.12  •  03 

Notes: 

Te  =  electron  temperature;  =  Te/300®K. 


tti:  Branching  ratios  from  Guberman  [1988]  for  i;"  =  0  state;  dayside  mean  vibrational  population  not  reliably 
known. 

0^3:  Branching  ratios  from  Guberman  [1991]  for  t;"  =  0  state;  dayside  mean  vibrational  population  not  reliably 
known.  See  Queffelec  ei  al  [1985]  for  laboratory  studies  indicating  a  total  N(^D®)  yield  ~  1.8.  See  Noren 
et  al  [1989]  for  recent  laboratory  results  suggesting  a  need  for  revision  of  the  adopted  rate. 
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Table  A3:  AURIC  Metastable- Electron  Quenching  Rates 


Quenching  Channel 

Rate  Coefficient  (cm^  s“^) 

References 

0+(2po)^g- 

0+('‘5®)  +  e- 

qia  =  3.31  X  10-®re"^^^ 

Pradkan  [1976] 

0+(2D‘’)  +  e- 

qn  =  1.39  X  10- 

0+(2po)^g- 

0+('‘5»)  +  e- 

32  =  6.60  X  10-®re"^^^ 

N(2p“)  +  e- 

N('‘5®)  +  e- 

33a  =  1-60  X  10-12  2^0.85 

Berrington  and  Burke  [1981] 

N(2Z)")  +  e- 

33j  =  9.50  X  10-® 

N(2D»)  +  e- 

N(‘‘5“)  +  e- 

34  =  3.80  X  10-12  T°  ®1 

0(i5)  +  e- 

0(3p)  +  e- 

35„  =  7.30  X  10-13  2;o.94 

Berringion  and  Burke  [1981] 

— +  0(  ^D)  +  e~ 

35i  =  8.50  X  10-3 

0(  i£»)  +  e- 

0(3p)  +  e- 

36  =  1.60  X  10-12  T®-®! 

Notes: 

=  electron  temperature;  Tq  =  Te/300®K. 

Here  and  in  Table  A4,  designation  for  ground  state  oxygen  ions  is  used  for  clarity;  elsewhere,  the 

designation  is  simply  O'*" . 
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Table  A4:  AURIC  Radiative  Decay  Channels 


Channel 

A  Coefficient  (s~^) 

References 

0+(2po) 

-*  0+(^5“)  +  H2470A) 

>12470  =  4.53  X  10-2 

Zeippen  [1982] 

0+(2£>®)  +  /ii/(7320A) 

A7320  =  9.85  X  10-2 

0+(2£)®)  +  /«/(7330A) 

>17330  =  7.49  X  10-2 

0+(2d<>) 

0+(‘‘5®)  +  fti/(3727A) 

2I3727  =  8.89  X  10-® 

N(2po) 

^  N('‘5®)  +  /ii/(3466A) 

.^3466  =  5.31  X  10-^ 

Butler  and  Zeippen  [1984] 

N(2D®)  +  /ij/(10400A) 

j4io4oo  =  8.05  X  10-2 

^  N(''5®)  +  /iz/(5200A) 

j452oo  =  1-28  X  10“® 

0(^5) 

0(^P)  +  /i2/(2972A) 

^2972  =  7.63  X  10-2 

Baluja  and  Zeippen  [1988] 

0(  ^D)  +  hv{f>bnk) 

-'45577  =  1.215 

0('D) 

0(3p)  +  /ii/(6300A) 

•'46300  ”  5.63  X  10  ^ 

0(3p)  +  fei/(6364A) 

'46364  =  1.82  X  10“^ 

N2(yl) 

N2(X)  +  /ii/(VK) 

^VK  =  0.422 

Piper  [1993b] 

Notes: 

Avk-  Assuming  N2(A)  in  the  v”  —  0  state;  dayside  mean  vibrational  population  not  reliably  known. 

Here  and  in  Table  A3,  designation  for  ground  state  oxygen  ions  is  used  for  clarity;  elsewhere,  the 

designation  is  simply  0+ . 
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ABSTRACT 

Computational  Physics,  Inc.  (CPI)  has  been  under  contract  to  the  Air  Force  Phillips  Laboratory/Geophysics 
Directorate  to  re-engineer  CPI  dayglow  and  nightglow  codes  conforming  to  modem  programming  standards, 
followed  by  testing  and  validation  of  these  codes.  The  software  and  associated  input  files  comprise  that  portion  of 
AURIC  for  the  modeling  of  dayglow  and  nightglow  from  100  to  about  900  nm  (MODTRAN  comprises  the  other 
portion).  This  paper  will  describe  the  current  capabilities  of  the  dayglow  portion  of  AURIC  and  present  selected 
dayglow  results.  The  dayglow  portion  models  photoionization,  photodissociation,  photoelectron  energy  degradation, 
photoelectron  impact  excitation,  chemistry,  solar  resonance  and  resonant  fluorescence  scattering,  and  multiple 
scattering  of  optically  thick  lines.  Calculated  observables  are  ion  and  neutral  densities  and  radiances  of  atomic  lines 
and  molecular  band  systems.  AURIC  is  currently  being  used  to  analyze  satellite  optical  data,  perform  studies  of 
production  and  loss  processes,  and  support  satellite  experiments  through  data  simulations. 

Keywords:  UV-visible  optical  backgounds,  remote  sensing,  dayglow 


1.  INTRODUCTION 

This  paper  discusses  the  model  AURIC  that  has  been  under  development  by  Computational  Physics,  Inc. 
(CPI)  for  the  Air  Force  Phillips  Laboratory/Geophysics  Directorate  (PL/GP).  A  more  detailed  version  of  the  paper' 
is  in  progress  for  submission  to  J.  Geophys.  Res,  That  paper  will  include  descriptions  of  both  dayglow  and  nightglow 
capabilities,  a  complete  reference  list,  a  more  detailed  mathematical  description  of  the  photoelectron  flux  model,  and 
details  about  electron  impact  cross  sections  and  rate  coefficients.  Under  the  work  being  documented,  an  AURIC 
User’s  Manual  has  also  been  written\  The  manual  contains  detailed  information  on  how  to  run  AURIC  on  Unix, 
VMS,  and  PC  systems  and  includes  numerous  tables  and  figures  of  inputs  and  outputs.  AURIC  stands  for 
Atmospheric  Ultraviolet  Radiance  Integrated  Code;  the  term  Integrated  refers  to  the  coupling  of  the  suite  of  models 
to  be  discussed  here  with  MODTRAN.  The  integrated  package  allows  one  to  calculate  spectral  radiances  that 
contain  dayglow  and  Rayleigh  scattering  of  sunlight  (where  appropriate)  or  nightglow  and  Rayleigh  scattering  of 
moonlight  (again,  where  appropriate).  MODTRAN’ s  overall  capabilities  include  rapid  molecular  band  model 
calculations  of  radiances  in  the  IR,  transmittances  from  the  IR  to  the  UV,  and  Rayleigh  scattering  of  sunlight  and 
moonlight^.  Early  work  by  CPI  was  reported  in  a  1992  SPIE  paper  by  Link  et  al.^  The  models  being  provided  by 
CPI  calculate  dayglow,  nightglow,  and  an  assortment  of  densities  above  90  km.  The  dayglow  portion  of  the  software 
has  its  origin  in  codes  developed  by  CPI  and  the  Naval  Research  Laboratory  (NRL)  under  the  name  PEGFAC*.  The 
PEGFAC  software  has  been  rewritten  using  modem  programming  standards.  In  addition,  major  improvements  have 
been  made  in  the  handling  of  I/O.  Beyond  these  improvements,  the  dayglow  portion  of  AURIC  also  has  many 
enhancements  over  PEGFAC.  Examples  include  the  addition  of  a  global  H  Lyman  a  model  (contributed  by  one  of 
us  -  JEB)  and  major  revisions  to  1)  various  sets  of  electron  impact  cross  sections,  2)  rate  coefficients,  and  3) 
synthetic  spectra  of  various  molecular  band  systems.  The  nightglow  portion  of  AURIC  also  has  its  origins  in  CPI 
software  that  has  been  rewritten  in  the  same  spirit  as  the  dayglow  software.  There  are  two  modules.  One  calculates 
emission  from  the  F2  region  of  the  ionosphere  (arising  from  recombination  of  O  and  Oj  with  electrons);  the  other 
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calculates  emission  from  the  upper  mesosphere/lower  thermosphere  (arising  from  O  atom  recombination).  The  total 
number  of  codes  calculating  daygiow  and  nightglow  is  27.  Included  is  a  code  that  provides  a  text-based  user 
interface  to  executing  the  daygiow  and  nightglow  codes  and  another  that  serves  as  an  interface  to  MODTRAN.  The 
CPI  software  within  AURIC  is  comprised  of  approximately  35,(K)0  lines  of  code. 

The  primary  sponsor  of  the  work  described  in  this  paper  is  the  DoD  Ballistic  Missile  Defense  Office  with 
technical  monitoring  through  PL/GP.  The  interest  on  the  part  of  the  sponsor  is  to  provide  a  tool  to  the  optical 
backgrounds  community  for  characterizing  naturally  occurring  atmospheric  emissions  from  the  far  ultraviolet  (PTJV) 
to  the  near  infrared  as  seen  from  space.  As  noted  above,  present  capability  allows  for  the  characterization  of 
emission  on  both  the  dayside  and  nightside.  Future  capability  will  include  the  modeling  of  auroral  and  twilight 
emissions.  In  addition  to  characterizing  daygiow  and  nightglow,  AURIC  can  also  be  used  to  perform  analyses  of 
satellite  and  rocket  data  and  studies  of  a  variety  of  neutral  and  ion  species  on  the  dayside.  An  example  of  such  a 
study  was  recently  published  by  Siskind  et  al.’  addressing  the  effects  of  enhanced  solar  soft  X-ray  fluxes  on  the 
density  of  nitric  oxide  (NO)  in  the  lower  thermosphere. 

Table  1  lists  atomic  daygiow  and  nightglow  features  presently  contained  in  AURIC.  Similar  information 
for  band  systems  is  given  in  Table  2.  Emission  for  some  of  the  daygiow  and  nightglow  features  arises  from 
metastable  states.  Densities  of  these  states  are  calculated  by  AURIC’s  photochemistry  and  nightglow  chemistry 
models.  Table  3  lists  the  daygiow  metastables  which  are  part  of  a  larger  list  that  includes  ground  state  ion  species 
and  N("S). 


Table  1.  Atomic  lines,  transitions,  type  of  airglow  (day/night)  modeled,  and  identification  of  those  features  affected 
by  scattering  or  requiring  chemistry  modeling  (metastables) 


Species 

Transition 

Wavelength 

(nm) 

Region 

Comment 

H 

p  s 

121.6 

Global 

Requires  multi-dimensional  photon 

transport  model 

N 

'P^  D 

149.3 

Day 

N 

2  .2  o 

p-»p 

174.3 

” 

N 

2  0.4  0 

p  s 

346.6 

Metastabie 

N 

2  0  .4  0 

D^S 

520.0 

Metastable 

N+ 

3  0.3 

D  P 

108.5 

Day 

N+  . 

S  0,3 

S-^P 

214.3 

It 

Metastabie 

0 

130.4 

Day/Night 

Requires  photon  transport  model 

0 

135.6 

Requires  photon  transport  model 

0 

S“^  P 

297.2 

Metastabie 

0 

'S^D 

557.7 

Metastabie 

0 

630.0 

Metastabie 

0 

777.4 

0 

'p-^V 

844.6 

0+ 

2  0  .4  o 

p-^S 

247.0 

Day 

Metastabie 

0+ 

2  0.4  0 

D  S 

372.7 

Metastabie 

0+ 

2  0,2  o 

P  D 

732.0 

Metastabie 
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Table  2.  Band  systems,  transitions,  type  of  airglow,  and  wavelength  regions  covered. 


Species 

Transition 

Band  system 

Region 

Wavelength  range  (nm) 

N2 

a  n<y  "4  X 

Lyman-Birge*Hopfield 

Day 

125  -  240 

N2 

3  +  1  + 

AZu  -^XZo 

Vegard-Kaplan 

150  -  690 

N2 

CDu  ->Bng 

Second  Positive 

•* 

280  -  460 

N2 

B  Og  -4  A  Zjj 

First  Positive 

600  to  beyond  1000 

N2+ 

2  +  2  + 
BZu  -»XZ„ 

First  Negative 

Day 

360  -  520 

N2+ 

A'Du  X  Zo" 

Meinel 

500  to  beyond  1000 

NO 

2^+  2  ® 

AZ  -»xn 

Y 

Day 

200  -  310 

NO 

c^n  ^  x^n 

2  +  2 

5 

ft 

190  -  250 

NO 

DZ  -»xn 

8 

180  -  260 

02 

Atmospheric 

Night 

680  -  1000 

02 

/3  I 

A  Ay  — >  a  Ag 

Chamberiain 

300  -  440 

02 

A^z;^  -4X^Ze' 

Herzberg  I 

" 

250  -  480 

02 

C  Zy  -4  X  Zo” 

Herzberg  II 

260-470 

Table  3.  Species  treated  by  the  photochemistry  module  and  identification  of  dominant  production  mechanisms 


Species 

Dominant  Production  Mechanisms 

Slow  reacting  species  (time  integrated) 

no" 

Charge  transfer 

02* 

photoionization,  charge  transfer 

0* 

e  impact,  photoionization 

Nf's) 

N(^)  quenching  by  O 

Fast  reacting  species 

n"" 

e  impact,  photoionization 

n; 

+  2 

e  impact,  photoionization 

0(D) 

2 

e  impact,  photoionization 

0(P) 

e  impact,  photoionization 

nCd) 

e  impact,  photodissociation 

N(P) 

1 

e  impact,  photodissociation 

0(D) 

J 

e  impact,  O2  recombination 

o(S) 

e  impact,  quenching  of  N2(A)  state 

N3(A) 

e  impact 
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2.  SOFTWARE  ivIODULES 


Figure  1  shows  a  flow  diagram  that  contains  the  AURIC  dayglow  modules.  Table  4  gives  a  brief 
description  of  each  module.  The  basic  steps  are:  (1)  specifying  needed  inputs  that  include  parameters  for  specifying 
a  model  atmosphere,  ionosphere,  and  solar  EUV  flux,  (2)  evaluating  these  latter  inputs,  (3)  calculating  photoelectron 
fluxes,  (4)  calculating  volume  production  rates  by  electron  impact,  photoionization,  photodissociation,  and  resonant 
fluorescence,  (5)  calculating  densities  affected  by  chemistry,  (6)  calculating  intensity  profiles  for  disk-to-limb 
viewing,  and  (7)  representing  these  intensities  in  spectral  radiance  form  (Rayleighs  nm  *)  using  molecular  synthetic 
spectral  routines. 


BEGIN 


Figure  1.  Flow  diagram  of  modules  contained  in  dayglow  portion  of  AURIC. 
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Table  4.  Brief  description  of  dayglow  modules  appearing  in  Figure  1. 


Module 

Description 

GEOPARM 

Calculates  selected  geophysical  parameters 

ATMOS 

Specifies  model  atmosphere  (MSISE-90  4-  SHARC 
atmospheric  generator) 

lONOS 

Specifies  model  ionosphere  using  FAIM 

SOLAR 

Specifies  solar  EUV  flux  using  Hinteregger  algorithm 

GOLDEN 

Calculates  slant  column  densities  of  N2,  0^,  0,  NO,  and 
O3  for  given  SZA 

RESOURCE 

Calculates  photoelectron  source  functions 

PEFLUX 

Calculates  photoelectron  fluxes 

e.IMPACT 

Calculates  volume  production  rates  due  to  electron  impact 
excitation 

DAYCHEM 

Calculates  ion  and  neutral  densities 

SUMVER 

Places  all  volume  emission  rates  in  a  single  file 

LOSDEN 

Calculates  line-of-sight  column  densities  of  and  NO 

needed  for  calculation  of  N.,*  ING,  N2^  Meinel,  and  NO 
y,  5,  £  band  intensities 

LOSINT 

Performs  line-of-sight  integrations  to  obtain  intensities 

LYMAN 

Calculates  HI  121.6  nm  intensity 

REDISTER 

Calculates  intensities  of  OI  130.4  nm  and  01  135.6  nm 
features 

Synthetic  spectral  routines  for  molecular 

Distributes  total  molecular  band  system  intensities  among 

band  systems 

appropriate  bands 

RAYSCAT 

Converts  MODTRAN  output  to  AURIC  synthetic  spectral 
file  format  for  Rayleigh  scattering  of  sunlight 

MERGESYN 

Merges  all  synthetic  spectra 

3.  DA YSroE  MODELING 


Non-resonant  excitation  mechanisms 
Solar  EUV  model 


The  states  associated  with  dayglow  and  species  introduced  in  Section  1  are  produced  by  one  or  more  of  the 
following  processes:  photoionization,  photodissociation,  solar  resonance  and  resonant  fluorescence  scattering, 
photoelectron  impact,  photochemistry,  and  multiple  scattering  for  optically  thick  lines.  The  starting  point  for  all  of 
these  processes  is  the  unattenuated  solar  flux  as  a  function  of  wavelength  X,  The  wavelength  range  of  interest  is 
from  -100  nm  down  to  X-ray  wavelengths.  The  shortest  wavelength  presently  treated  by  AURIC  is  1  nm.  Most  of 
the  solar  energy  deposited  at  shorter  wavelengths  occurs  below  100  km  and  is  insignificant  except  under  extreme 
solar  storm  conditions.  AURIC  uses  a  discrete  version  of  the  solar  spectrum  from  Hinteregger  et  al.‘®  extended  by  us 
from  1.8  down  to  1  nm  (see  Siskind  et  al.’). 
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Photoelectron  model 


AURIC  calculates  the  photoelectron  flux  in  the  local  approximation.  Numerous  local  models  have  been 
published  over  the  past  20  years.  Examples  are  Victor  et  al."  ,  Jasperse’%  Richards  and  Torr‘\  and  Strickland  and 
Meier*.  The  solution  is  adequate  for  most  applications  involving  photochemistry  and  dayglow.  Transport  effects 
should  be  taken  into  account  if  one  is  concerned  with  dayglow  above  about  300  km.  The  usual  situation  where  this 
concern  exists  is  in  modeling  twilight  photochemistry  and  emissions.  The  form  of  the  local  equation  treated  by 
AURIC  is: 


0  =  “X  n;(z)Of  (£’)(t)(z,  E) 

t 

+  E)(|)(Z,E  ’) 

i  ^ 

+  n,(z)^[4(z,E)(|)(z,E)] 


(l) 


Terms  include: 

<j)  spherical  electron  flux  (e  cm  V^eV**) 

n£  density  of  species  (N^,  0„  or  O) 

total  cross  section  of  species  (ionization  plus  excitation) 

R^  energy  redistribution  function  (eV  ) 

Lg  loss  function  for  energy  loss  to  thermal  electrons. 

source  function  for  photoelectron  production  by  solar  EUV  flux  (e  cm  s  eV  ). 

The  dependence  on  the  solar  zenith  angle  (SZA)  is  not  displayed  in  eq.  1  nor  in  the  remaining  equations  in  this 
section  in  order  to  simplify  their  appearance.  The  term  explicitly  dependent  on  SZA  is  the  source  function.  This 
function  contains  contributions  from  both  photoelectrons  and  Auger  electrons.  The  form  of  R^  is: 


R,(E  E)  =-^ - — - .  eV-l  (2) 

OfiE) 

where  the  sum  is  over  ionization  and  excitation  processes  and  contains  differential  cross  sections  in  cm  eV  .  Eq.  1 
is  solved  on  an  electron  energy  grid  from  1  to  800  eV.  Strickland  and  Meier  discuss  the  numerical  technique  used 
to  solve  the  equation. 

Volume  production  rates 

Here,  we  address  production  rates  by  the  following  processes:  photoionization,  photodissociation,  and 
electron  impact  Rates  for  the  first  two  processes  are  obtained  from: 
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-3  -1 

cm  s 


(3) 


Pk.ik  (2)  =  «f  (z)a,^*""’  (A, ) 


where  tcF  is  the  attenuated  solar  EUV  flux  at  altitude  z  (accounting  for  slant  paths  associated  with  a  non-vertical  sun) 
and  c  is  the  photo-cross  section  for  process  k  involving  species  1.  A  similar  expression,  generalized  to 
integrations  over  individual  lines,  is  used  for  resonant  fluorescence  scattering.  The  latter  process  is  of  interest  within 
the  AURIC  model  for  calculating  scattering  within  the  bands  of  the  NO  y,  6,  and  e  systems  as  well  as  the  ING 
and  Meinel  systems.  Eq.  3  is  used  to  calculate  production  rates  for  the  species  in  Table  3  and  the  photo-component 
of  the  Nn  108.5  nm  emission  rate. 

Production  by  electron  impact  is  given  by: 


^rii; 

J 


-3  -1 
cm  s 


(4) 


Many  of  the  emission  features  within  the  AURIC  model  are  produced  exclusively  by  electron  impact.  Production 
rates  needed  for  chemistry  modeling  are  typically  comprised  of  terms  from  both  eqs.  3  and  4. 

Photochemistry 

Table  3  lists  the  species  for  which  densities  are  calculated  by  AURIC’ s  photochemistry  model.  The 
densities  of  the  species  labeled  as  slow  reacting  are  calculated  by  solving  coupled  rate  equations  of  the  following 
form: 


tin- 

=  Pi(Zyt)  -  li(z.t)  cmV  (5) 

where  p,  and  L  are  the  total  production  and  loss  rates.  The  derivative  is  set  to  zero  for  calculating  the  densities  of  the 
fast  reacting  species.  The  possible  terms  comprising  a  given  total  production  rate  are  rates  from  eqs.  3  and  4  and  one 
or  more  rates  involving  chemical  reactions.  Loss  processes  include  quenching,  spontaneous  emission,  and  reactions 
such  as  atom-ion  interchange.  A  detailed  description  of  the  chemical  scheme  with  a  complete  listing  of  reactions  and 
associated  rates  is  given  in  reference  1 . 


Solar  resonance  and  resonant  fluorescence  scattering  features 

Features  to  be  addressed  are  01 130.4  nm,  01  135.6  nm,  HI  121.6  nm,  the  y,  5,  and  £  systems  of  NO,  and 
the  first  negative  and  Meinel  systems  of  N2^.  The  01  130.4  nm  "feature"  is  actually  a  triplet  (130.2  nm,  130.5  nm, 
130.6  nm)  arising  from  the  transition  (2p)  ^2.1.0  ^  The  initial  excitation  rate  of  the  state,  denoted 

S^,  includes  both  electron  impact  and  solar  resonant  absorption  terms.  For  the  optically  thick  01 130.4  nm  triplet,  as 
well  as  for^the  much  more  optically  thin  OI  135.6  nm  feature,  the  transport  equation  formulation  and  solution  code  of 
Gladstone  is  used.  This  code,  named  REDISTER,  employs  the  Feautrier  technique  to  solve  the  radiative  transfer 
equation  for  resonance  multiplet  emissions  in  a  plane-parallel  planetary  atmosphere.  Frequency  redistribution  is 
handled  in  the  angle-averaged  partial  redistribution  approximation,  taking  full  account  of  vertical  variations  in 
temperature,  scatterer  density,  and  the  densities  of  absorbing  species  (as  provided  by  the  MSISE-90‘'  atmosphere 
model).  The  initial  internal  (i.e.,  photoelectron  impact)  source  term  is  provided  by  the  e_IMPACT  module,  while  the 
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external  (i.e.,  solar  resonant  absorption)  source  term  is  evaluated  within  REDISTER.  Output  comprises  optical 
depths,  total  source  functions  S,  and  line-of-sight  intensities. 

The  Lyman  a  (HI  121.6  nm)  radiance  calculations  within  AURIC  are  performed  by  the  LYMAN  module, 
which  is  based  on  a  nonisothermal  adaptation  of  the  complete  frequency  redistribution  algorithm  of  Anderson  and 
Hord‘*.  As  implemented  in  AURIC,  it  takes  into  account  the  variation  of  densities  and  temperature  with  altitude  but 
assumes  spherically  symmetric  distributions.  The  temperature,  H  density,  and  (the  absorbing  species)  density 
profiles  are  taken  from  MSISE-90*^  for  altitudes  74  -  455  km;  at  higher  altitudes,  an  analytic  formulation’’  is  used  to 
estimate  H-atom  densities.  The  default  incident  solar  line  center  Lyman  a  flux  is  estimated  using  the  same 
scaling  of  the  reference  Hinteregger  SC21REFW  spectrum  as  is  used  for  specifying  the  incident  solar  EUV  flux, 
taking  the  integrated  and  line  center  fluxes  to  be  numerically  equal.  The  anisotropic  scattering  phase  function  is 
included  in  the  calculation  of  the  single  scattering  component  of  the  total  intensity. 

Resonant  fluorescence  g-factors  for  the  y,  5,  and  e  band  systems  of  NO  are  currently  computed  using  scaled 
solar  fluxes  from  Mount  and  Rottman’*  with  oscillator  strengths  from  Huber  and  Herzberg’’,  as  described  in  Cleary’®. 
NO  densities  are  obtained  from  an  empirical  model’’.  The  total  system  g-factors  for  the  y  (v’  =  0  -  4  progressions),  5 
(v’  =  0  progression),  and  £  (v’  =  0  -  2  progressions)  are  2.9  x  10  s  ,  3.1  x  10  s  ,  and  5.8  x  10  s  ,  respectively. 
For  the  ING  (first  negative)  system,  optically  thin  g-factors  have  been  computed  using  the  1985  high  resolution 
SUSIM  solar  irradiance  spectrum"  and  the  molecular  constants  compiled  by  Laux  and  Kruger‘S*  the  total  system  g- 
factor  is  0.148  s  .  The  g-factor  for  the  N^^  Meinel  system  adopted  for  AURIC  is  0.114  s  ,  derived  using  the 
oscillator  strengths  of  Cartwright‘S  and  the  MODTRAN  reference  solar  irradiance  spectrum.  Opacity  and  albedo 
effects  are  currently  ignored  for  all  these  systems. 


Dayglow  spectral  radiances 
The  intensity  equation 

AURIC  calculates  spectral  radiances  in  units  of  Rayleighs  nm  *  for  user  specified  look  angles  from  0^ 
(zenith)  to  180*^  (nadir).  The  emission  is  comprised  of  molecular  bands  and  atomic  lines.  Rotational  lines  within 
these  bands  along  with  the  atomic  lines  are  smoothed  to  the  user  specified  wavelength  resolution.  The  spectral 
radiance,  designated  by  47clj^  is: 

4k /x  =  W^J  J(s,X)T(SyX)e''^^^'^^ds  Rayleighs  nm’  (6) 

with  terms  as  follows: 

j  volume  emission  rate  in  photons  cm  ^s  *nm  ^ 

T  transmission  function  for  self  absorption  (dimensionless) 

t  pure  absorption  optical  depth  (dimensionless) 
s  distance  from  observing  point  along  line-of-sight. 

Line-of-sight  integration  is  through  a  spherically  stratified  atmosphere.  The  solar  zenith  angle  is  assumed  constant 
along  the  line-of-sight  and  thus,  the  present  version  of  AURIC  is  not  suitable  for  calculating  intensities  near  the 
terminator. 

The  transmission  function  T  is  unity  for  most  of  AURIC’s  emission  features.  Exceptions  are  01  130.4  nm, 
01  135.6  nm,  and  HI  121.6  nm.  The  T  functions  for  these  features  are  treated  by  the  transport  models  described  in 
the  previous  section.  The  exponential  term  is  unity  long  ward  of  -175  nm.  Below  this  wavelength,  t  is  non-zero  due 
to  Schumann-Runge  continuum  absorption  by  O2  which  peaks  at  -140  nm. 
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Spectral  radiances  in  nm  are  not  calculated  by  AURIC  from  eq.  6.  Instead,  total  intensities  by  line,  band, 
or  band  system  are  first  calculated  followed  by  separate  procedures  for  distributing  these  intensities  in  wavelength 
for  the  selected  resolution.  More  details  follow  below  about  the  synthetic  spectral  routines  that  produce  the 
wavelength  dependence  of  bands  and  band  systems.  From  eq.  6,  the  total  intensity  of  a  line,  band,  or  band  system  is: 


471  /,  =  10 


Rayleighs 


(7) 


where  T  has  been  set  to  unity  (applicable  except  for  those  features  treated  by  the  transport  models),  the  subscript  i 
refers  to  a  line,  band,  or  system,  and  units  of  j.  are  photons  cm  s  \  Three  forms  of  j  are  treated  by  AURIC.  The  first 
is  given  by  eq.  4  for  electron  impact  excitation  with  unit  probability  of  emission.  The  second  is  given  by: 


ji(z)  =  A,n(z) 


(8) 


which  is  appropriate  to  metastable  states  where  n  is  the  density  of  a  given  state  calculated  by  AURIC’ s 
photochemistry  model  and  A^  is  the  transition  rate  associated  with  the  given  feature.  The  third  form  refers  to  solar 
resonant  fluorescence: 


Ji  =  gMz) 


(9) 


where  g.  is  the  g-factor  discussed  in  the  previous  section.  This  form  of  j  is  used  to  calculate  the  total  volume 
emission  rates  for  the  NO  and  systems. 

Molecular  band  synthetic  spectral  models 

To  incorporate  molecular  emissions  in  total  atmospheric  radiance  spectra,  it  is  necessary  to  distribute  the 
system  emission  rate  among  the  observable  bands;  in  turn,  the  distribution  of  emission  within  a  single  band  is 
dictated  by  the  rotational  line  spectrum.  For  the  band  systems  listed  in  Table  2,  spectral  synthesis  routines  have  been 
adopted  from  numerous  sources.  In  all  cases,  spectra  are  evaluated  at  0.1  nm  resolution  and  subsequently  smoothed 
to  the  user-specified  resolution.  Spectral  synthesis  routines  for  the  bands  of  the  NO  y,  5,  and  e  systems  are  from 
Cleary" ,  with  excitation  out  of  the  ground  state  v"  =  0  level.  For  the  N^  band  systems,  spectral  modeling  of 
individual  bands  is  earned  out  using  routines  provided  by  D.  Gattinger.  Individual  band  intensities  have  been 
evaluated  using  the  solar  flux  spectra  and  transition  probabilities  noted  previously,  coupled  with  an  assumed  N,^ 
ground  state  vibrational  temperature  of  6000  K  which  gives  a  relative  population  of  0.40  for  the  v"=  0  level. 

Modeling  of  the  photoelectron  impact-produced  N^  LBH,  VK,  and  2PG  emission  spectra  utilizes  the 
spectral  modeling  routines  of  Conway^  and  Conway  and  Christensen"^  For  the  LBH  system,  progressions  v’  =  0  -  6 
are  treated  with  the  relative  populations  estimated  from  the  v -0  Franck-Condon  factors,  taking  into  account 
wavelength-dependent  photoabsorption  by  Progressions  v’  =  0  -  9  are  included  in  the  modeling  of  the  VK 
system,  using  the  relative  populations  of  Cartwright^.  This  system  is  affected  by  quenching;  in  AURIC,  the 
dependence  of  the  quenching  rate  on  v'  is  currently  ignored.  Modeling  of  the  2PG  system  includes  progressions  v’  = 
0  -  4,  again  with  the  v'  relative  populations  estimated  from  the  v'-O  Franck-Condon  factors.  The  IPG  emission 
system  is  also  generated  by  photoelectron  impact.  For  this  system,  spectral  modeling  of  individual  bands  is  done 
with  a  routine  provided  by  D.  Gattinger;  the  relative  distribution  of  the  total  system  intensity  among  system  bands  is 
taken  from  Table  4.12  of  Vallance  Jones“. 
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Rayleigh  scattering  of  sunlight 

An  interface  to  MODTRAN  exists  for  calculating  Rayleigh  scattered  spectral  radiances,  saving  them  to  a 
file,  reading  the  file  at  an  appropriate  point  within  the  execution  of  the  sequence  of  AURIC  modules,  and  adding 
these  radiances  to  dayglo w  spectral  radiances.  Rayleigh  scattering  must  be  taken  into  account  longward  of -180  nm 
at  tangent  altitudes  below  -100  km.  MODTRAN’s  full  capability  includes  rapid  molecular  band  model  calculations 
of  radiances  in  the  IR,  tranmittances  from  the  IR  to  the  UV,  and  Rayleigh  scattering  of  sunlight  and  moonlight  (see 
MODTRAN  references  in  Section  1).  The  user  must  provide  his/her  own  copy  of  MODTRAN  for  interfacing  with 
that  portion  of  AURIC  being  discussed  in  this  paper. 

Key  inputs  for  dayside  modeling 

Table  5  lists  the  key  inputs  to  AURIC’s  dayglow  modules.  The  first  five  entries  are  obtained  from  models 
or  algorithms  as  identified  in  the  Reference  column.  The  remaining  entries  refer  to  a  large  volume  of  data  obtained 
from  compilations  by  either  the  authors  or  other  investigators.  To  conserve  space,  we  refer  the  reader  to  the  papers 
referenced  in  the  table.  Most  of  the  data  appear  explicitly  in  several  tables  appearing  in  Strickland  et  al*  (for 
example,  extensive  tables  of  rate  coefficients  for  AURIC’ s  dayside  and  nightside  chemistry  models). 


Table  5.  Key  inputs  to  dayglow  modules  and  references. 


Parameter 

Reference 

Model  atmosphere  (N,,  0,,  and  0) 

MSISE-90  (Hedin“) 

Model  atmosphere  (NO) 

Smith  et  al.'‘ 

Model  ionosphere 

FAIM  (Anderson  et  al.^) 

Geomagnetic  field  (for  relating  geomagnetic  to 
geographic  coordinates) 

IGRF(Peddie*) 

Solar  EUV  spectrum 

Hinteregger  et  al.‘® 

Photoabsorption  cross  sections 

Conway” 

Photoionizaiion  cross  sections 

Conway” 

Electron  impact  energy  loss  cross  sections 

Majeed  and  Strickiand^^ 

Electron  impact  cross  sections  for  calculating 
production  rates  needed  by  chemistry  model 

Strickland  et  al.‘ 

Electron  impact  cross  sections  for  calculating 
emission  rates 

Strickland  et  al.‘ 

Schumann-Runge  0,  continuum  photo¬ 
absorption  cross  section 

Hudson” 

Rate  coefficients 

Strickland  et  al.‘ 

Transition  rates 

Strickland  et  al.‘ 

Population  distributions 

Strickland  et  al.‘ 
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4.  SELECTED  DAYGLOW  RESULTS 


^  A  detailed  description  of  the  output  from  AURIC’ s  day  glow  modules  is  available  in  the  AURIC  User’s 
Manual'.  A  subset  of  some  of  the  more  interesting  results  will  be  illustrated  here.  Results  will  be  shown  in  the  form 
of  photoelectron  fluxes,  altitude  profiles  of  volume  emission  rates,  limb  intensity  profiles,  and  spectral  radiances  for 
selected  tangent  altitudes  as  seen  from  satellite  altitudes.  Calculations  have  been  performed  for  an  MSISE-90 
midlatitude  model  at  moderate  solar  activity  =  150)  for  a  solar  zenith  angle  of  60^  The  exception  to  this  was  the 
calculation  of  the  photoelectron  fluxes  for  which  conditions  were  chosen  to  match  those  of  the  data  to  which  the 
fluxes  will  be  compared.  Figure  2  shows  calculated  and  measured  photoelectron  fluxes  at  an  altitude  of  167  km. 
The  data  were  taken  from  the  paper  of  Lee  et  al."*  and  have  been  divided  by  1.5  to  reproduce  the  overall  magnitude 
of  the  model  results.  The  accuracy  of  the  data  has  been  discussed  in  various  papers  over  the  past  several  years. 
Richards  and  Torr*^  argue  that  the  fluxes  have  the  correct  magnitude  while  Strickland  and  colleagues  argue  that  the 
fluxes  should  be  reduced  by  a  factor  between  1.5  and  2.0  (e.g.,  Strickland  and  Anderson^O-  Conway^  presents 
independent  information  that  also  argues  for  a  reduction  in  the  measured  photoelectron  fluxes.  The  different  points 
of  view  arise  from  the  use  of  different  sets  of  N,  excitation  cross  sections.  More  details  are  given  in  the  paper  by 
Majeed  and  Strickland^l  Returning  to  Figure  2,  two  model  spectra  are  shown  for  different  solar  EUV  fluxes  below 
250  A.  A  better  fit  to  the  data  is  obtained  by  increasing  the  Hinteregger  fluxes  at  short  wavelengths.  We  have 
chosen  a  scaling  factor  of  2  to  demonstrate  this  point.  Richards  and  Torr’^  have  also  argued  for  higher  fluxes  at 
shorter  wavelengths  based  on  the  same  types  of  shape  comparisons  between  their  model  results  and  the  Lee  et  al. 
data. 


Figure  2.  Calculated  and  measured^  photoelectron  fluxes  at  167  km  (the  latter  scaled  downward  by  a  factor  1/1.5). 
The  solid  curve  shows  the  effect  of  increasing  the  solar  flux  below  25  nm  by  a  factor  of  two. 
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Fi^  jre  3  illustrates  the  altitude  dependence  of  volume  emission  rates  for  selected  key  emission  features.  For 
illustration,  we  have  selected  N,  LBH,  IPG,  N,  VK  (includes  the  effect  of  quenching),  and  the  lines  of  atomic 
oxygen  at  135.6  nm  and  844.6  nm.  The  secondary  peak  at  low  altitudes  is  due  to  the  soft  X«ray  portion  of  the  solar 
spectrum.  Figure  4  shows  limb  profiles  corresponding  to  the  results  in  Figure  3  plus  profiles  for  the  optically  thick 
features  OI  130.4  nm  and  HI  121.6  nm.  As  noted  earlier,  the  REDISTER  model  of  Gladstone'^  is  used  to  calculate 
130.4  radiances  (as  well  as  OI  135.6  nm)  while  a  global  CFR  model  by  one  of  us  (JEB)  is  used  to  calculate  121.6  nm 
radiances.  For  LBH,  two  profiles  are  shown  for  emission  between  140  -  149  nm  and  170  -  180  nm;  the  difference  in 
shape  is  caused  by  differences  in  0,  pure  absorption  between  the  selected  intervals.  Spectral  radiances  from  100  to 
1000  nm  at  tangent  altitudes  of  80,  lOO.  and  120  km  are  shown  in  Figure  5.  Rayleigh  scattering  results  are  included 
which  begin  to  exceed  the  dayglow  results  as  the  tangent  altitude  decreases  below  100  km.  Tables  1  and  2  list  the 
dayglow  features  that  have  contributed  to  the  spectra  in  Figure  5. 


1  10  100  1000 
Volume  Emission  Rate  (cm’^  s"') 


Figure  3.  Selected  volume  emission  rates  for  a  60"*  solar  zenith  angle  and  a  Hinteregger  solar  EUV  spectrum  based 
on  F,^,,  =  150.  The  135.6  nm  rate  is  the  initial  rate  referred  to  as  S^  in  the  above  section  under  the  heading 
Solar  resonance  and  resonant  fluorescence  scattering  features.  X-rays  are  responsible  for  the  secondary 
low  altitude  peaks. 
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Figure  4.  Limb  profiles  in  Rayleighs  corresponding  to  the  volume  emission  rates  in  Figure  3  plus  profiles  for  the 
optically  thick  01  130.4  nm  and  HI  121.6  nm  features.  The  differences  in  the  shapes  of  the  two  LBH 
profiles  are  due  to  differences  in  pure  absorption  by  between  these  bands. 


Figure  5.  SpecQ-al  radiances  at  a  wavelength  resolution  of  1.0  nm  for  tangent  altitudes  of  80,  100,  and  120  km. 
Rayleigh  scattering  of  sunlight  dominates  the  spectrum  at  80  km  longward  of  1 80  nm. 
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5.  DISCUSSION  AND  SUMMARY 


As  of  this  writing,  final  testing  and  validation  are  underway  of  the  software  modules  introduced  in  Section 
2.  An  example  of  validation  efforts  is  presented  in  Figure  2  where  AURIC  photoelectron  fluxes  are  compared  to  the 
data  of  Lee  et  al.^.  Throughout  the  project,  CPI  has  made  the  airgiow  portion  of  AURIC  available  to  various  DoD 
affiliated  organizations.  Executables  along  with  all  needed  input  data  files  and  selected  graphics  routines  have  been 
installed  on  Unix,  VMS,  and  PC  systems  as  new  versions  of  the  software  have  become  available.  The  AURIC  User’s 
Manual  has  continued  to  evolve  from  version  to  version  of  the  software  and  has  been  included  with  the  above 
installations.  Examples  of  organizations  that  have  been  using  AURIC  are  the  sponsor  (PL/GP),  the  Naval  Research 
Laboratory,  the  Applied  Physics  Laboratory,  and  selected  government  contractors. 

End-of-project  documentation  has  already  been  cited  in  earlier  sections.  The  key  documents  are  a  science 
manual  and  the  user’s  manual.  The  science  manual  is  near  completion  and  is  being  written  as  a  paper'  for 
submission  to  an  appropriate  journal.  This  SPIE  paper  is  an  abbreviated  version  of  this  paper.  Present  plans  are  to 
submit  the  full  paper  to  J.  Geophys.  Res.  within  the  next  two  months.  Another  important  documentation  effort  has 
.been  directed  to  the  large  number  of  energy  loss  cross  sections  used  to  calculate  photoelectron  fluxes.  A  paper  has 
been  completed  on  this  subject  and  submitted  to  J.  Phys.  Chem.  Ref.  Data  by  two  of  us  (TM  and  DJS)^\  The  paper 
was  written  to  document  our  sets  of  energy  loss  cross  sections  for  N„  0„  and  O  which  contain  many  cross  sections 
from  recent  energy  loss  measurements.  Overall  assessments  of  these  sets  are  made  through  comparisons  with  total 
cross  sections  obtained  from  transmission  measurements  and  comparisons  of  cross-section-based  loss  functions  with 
the  Bethe  formula  above  100  eV.  Discussion  includes  energy  loss  cross  sections  used  by  other  modelers  in  the  field. 

The  airgiow  portion  of  AURIC  should  prove  to  be  a  valuable  tool  to  investigators  interested  in 
characterizing  optical  backgrounds  at  rocket  and  satellite  altitudes,  simulating  and  analyzing  optical  data,  and 
studying  chemically  active  species  (ions  and  neutrals)  in  the  thermosphere.  The  suite  of  models  is  easy  to  use 
through  either  the  supplied  interface  or  command  files.  All  input  files  are  formatted  which  allows  for  easy 
modification  of  parameters  such  as  and  arrays  containing  information  such  as  input  densities.  Since  the  modules 
are  first-principles  rather  than  empirically  based,  the  user  can  conduct  detailed  investigations  of  the  interrelationships 
between  physical  quantities  such  as  the  solar  EUV  flux  and  composition  on  the  one  hand  and  emissions  and 
chemically  active  densities  on  the  other. 
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NEW  SURVEY  OF  ELECTRON  IMPACT  CROSS  SECTIONS  FOR 
PHOTOELECTRON  AND  AURORAL  ELECTRON  ENERGY  LOSS 

CALCULATIONS 

Tariq  Majeed  and  Douglas  J.  Strickland 

Computational  Physics,  Inc.,  2750  Prosperity  Avenue,  Suite  600,  Fairfax,  VA  22031 

Newly  surveyed  sets  of  energy  loss  cross  sections  are  presented  for  N2,  Oj,  and  O.  The 
work  was  motivated  by  a  number  of  new  electron  energy  loss  measurements  in  the  late 
19805  and  early  19905  and  recent  selected  review  articles.  Each  set  includes  a  total 
ionization  cross  section  and  excitation  cross  sections  that  correspond  to  all  important  non¬ 
ionizing  energy  loss  channels  for  that  species.  A  total  cross  section  for  each  species  is 
constructed  by  summing  the  elastic  scattering  cross  section  with  the  ionization  and 
excitation  cross  sections.  The  sum  is  compared  to  a  measured  total  cross  section  obtained 
from  electron  transmission  experiments.  Good  agreement  is  achieved  for  each  of  the 
three  species.  A  loss  function  is  also  constructed  for  each  species  and  compared  with  the 
Bethe  formula  above  100  eV.  Good  agreement  is  also  achieved  in  energy  loss  which  is 
dominated  by  ion  and  secondary  electron  production.  Fluxes  of  photoelectrons  and 
auroral  electrons  have  been  calculated  for  the  new  sets  of  energy  loss  cross  sections  as 
well  as  our  previous  sets.  No  substantial  differences  occur  using  the  new  description  of 
energy  loss. 

Key  Words:  Energy  loss  cross  sections  for  N2,  O2  and  O,  Photoelectron  flux.  Auroral 
electron  flux,  Electron  impact.  Elastic  and  total  scattering  cross  sections.  Ionization  cross 
section,  Bethe  formula.  Inelastic  cross  section. 
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1.  INTRODUCTION 


The  work  described  below  was  undertaken  as  one  of  several  tasks  to  develop  a 
dayglow/nightglow  UV  radiance  model  for  the  integrated  model  AURIC  (Atmospheric 
Ultraviolet  Radiance  Integrated  Code).  The  term  integrated  refers  to  the  joining  of  this 
radiance  model  with  the  Air  Force  model  MODTRAN.  This  latter  model  provides  rapid 
molecular  band  model  calculations  of  radiances  in  the  IR,  transmittances  from  the  IR  to 
the  UV  and  Rayleigh  scattering  of  sunlight  and  moonlight^’^’^’".  The  designator  AURIC- 
R  will  be  used  to  distinguish  the  UV  radiance  portion  from  the  integrated  model. 
AURIC-R  has  been  under  development  by  Computational  Physics,  Inc.  (CPI)  over  the 
past  three  years  for  the  Geophysics  Directorate  of  the  Air  Force  Phillips  Laboratory 
(PL/GP).  A  key  task  has  been  the  re-engineering  of  Fortran  codes  within  the  PEGFAC 
(photoelectron  g-factor)  model^  using  modem  programming  standards.  Another  key  task 
has  been  I/O  restmcturing  and  updating  of  key  input  parameters.  Much  of  the  latter  effort 
has  been  directed  to  three  sets  of  electron  impact  cross  sections.  These  sets  are  used  to  1) 
perform  photoelectron  energy  loss  calculations,  2)  calculate  volume  production  rates  for 
chemistry  modeling,  and  3)  calculate  volume  emission  rates  for  specifying  spectral 
radiances.  This  paper  addresses  the  first  of  these  sets  containing  energy  loss  cross 
sections  for  N2,  O2,  and  O. 

The  motivation  for  this  work  comes  from  a  number  of  new  cross  section 
measurements  in  the  late  1980^  and  early  19905  (references  to  many  of  these 
measurements  appear  in  the  recent  reviews  of  Itikawa  et  al.^’^,  Itikawa  and  Ichimura*, 
Laher  and  Gilmore^  and  Kanik  et  al.  ^°).  Our  approach  has  been  to  gather  cross  sections 


41 


for  the  important  loss  channels  of  each  of  the  three  species  (a  single  channel  for 
ionication  and  several  for  excitation)  and  examine  them  in  two  ways.  First,  a  total  cross 
section  is  constructed  for  each  species  by  summing  the  total  ionization  and  energy  loss 
cross  sections  with  the  elastic  scattering  cross  section  for  that  species.  This  total  cross 
section  is  then  compared  with  measurements  from  electron  transmission  experiments. 
Second,  a  loss  function  is  constructed  and  compared  with  the  Bethe  formula  (see 
Strickland  et  al}^  for  its  form  and  application  to  Nj)  above  100  eV.  While  such  a  test  is 
not  useful  for  accurately  assessing  a  cross  section  set  at  energies  most  important  to 
photoelectron  energy  loss  calculations  (below  100  eV),  it  does  place  constraints  on  the 
total  inelastic  cross  section  and  differential  dependence  of  the  ionization  cross  section 
above  100  eV.  Assuming  good  knowledge  of  the  total  ionization  cross  section,  the 
differential  dependence  dictates  the  magnitudes  of  the  secondary  electron  energy  loss 
component  of  the  total  loss  function.  This  component  dominates  above  a  few  hundred  eV 
as  will  be  illustrated  later  in  the  paper. 

Ionization,  elastic,  and  total  cross  sections  of  Nj,  Oj,  and  O  appear  to  be  well 
quantified  at  this  time  through  both  laboratory  measurements  and  calculations.  As  we 
shall  demonstrate,  work  still  remains  to  be  done  in  quantifying  the  many  excitation  (non¬ 
ionizing  energy  loss)  channels  of  these  species,  especially  near  and  above  the  first 
ionization  threshold  (Rydberg  channels).  There  is  generally  good  agreement  among  the 
various  sets  of  total  ionization  cross  section  measurements  of  Nj,  O2,  and  O  (see  Itikawa 
et  al.^’^,  Itikawa  and  Ichimura®,  and  Kanik  et  al}^  for  specific  references).  The  most 
recent  measurements  of  total  ionization  cross  sections  of  N2  (Krishnakumar  and 


Srivastava'^)  and  O2  (Krishnakumar  and  Srivastava^^)  are  in  close  agreement  with  earlier 
measurements  of  Rapp  and  Englander-Golden*'^.  For  O,  the  recent  work  by  Itikawa  and 
Ichimura*  lends  support  to  the  measurements  of  Brook  et  al}^  The  total  and  elastic 
scattering  cross  sections  of  N2  and  O2  appear  to  be  well  characterized  based  on  agreement 
among  the  various  existing  sets  of  measurements  (see  reviews  of  Itikawa  et  a\^'^  and 
Kanik  et  for  agreement  among  the  original  data  sets).  For  O,  the  only  measurement 
of  the  total  cross  section  is  by  Sunshine  et  al}^  The  measurements  were  limited  to 
energies  from  1  to  100  eV  and  possess  more  scatter  than  for  the  corresponding  cross 
sections  of  N2  and  O2.  Itikawa  and  Ichimura®  constructed  a  total  O  cross  section  by 
summing  available  components  and  found  agreement  with  Sunshine  et  al}^  within  the 
scatter  of  the  data.  The  Itikawa  and  Ichimura®  cross  section  spans  a  larger  energy  range 
going  to  7000  eV.  Within  their  sum  is  an  elastic  scattering  cross  section  based  on 
calculations  rather  than  measurements.  The  available  measurements  are  by  Dehmel  et 
aiy’  which  appear  to  be  contaminated  by  inelastic  scattering®. 

Fox  and  Victor  discuss  electron  energy  loss  in  N2.  Cross  section  information  is 
in  the  form  of  loss  function  components  (for  excitation,  production  of  ion  states,  and 
kinetic  energy  of  secondaries)  with  direct  cross  section  information  limited  to  references. 
Several  compiled  sets  of  energy  loss  cross  sections  or  totals  by  species  have  been 
published  over  the  years  in  papers  addressing  the  calculation  of  photoelectron  and  auroral 
electron  fluxes  (e.g.,  Strickland  et  al}^,  Victor  et  al}^,  Oran  and  Strickland^®,  Jackman 
and  Green^*,  Mantas^^,  Stamnes  and  Rees^®,  Richards  and  Torr^'*’^^,  Solomon^®,  and 
Strickland  et  al.  ).  The  sources  of  measured  and  calculated  cross  sections  from  one  set 
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to  another  are  not  the  same  and  in  turn  can  lead  to  different  conclusions  from  analyses  of 
photoefjstron  data,  auroral  electron  data,  and  optical  data  involving  emission  features 
produced  by  electron  impact  excitation.  A  further  discussion  on  this  topic  will  be  given 
in  Section.  5. 

Since  this  paper  addresses  energy  loss  cross  sections,  there  will  be  limited 
discussion  of  the  collision  products  associated  with  a  given  loss  channel.  Collision 
products  are  important  for  the  other  two  cross  section  sets  mentioned  above,  namely  for 
production  rates  needed  in  chemistry  calculations  and  emission  rates  needed  in  radiance 
calculations.  Energy  loss  cross  sections,  on  the  other  hand,  are  used  to  calculate 
photoelectron  and  auroral  electron  fluxes  for  which  the  only  requirement  is  that  a  proper 
distribution  of  energy  loss  be  achieved  per  collision.  Here,  the  important  features  are 
excitation  thresholds,  cross  section  magnitudes,  and  in  the  case  of  ionization,  the  initial 
distribution  of  secondary  electrons  (for  a  discussion  of  the  treatment  of  secondaries  in 
AURIC-R  as  well  as  CPI’s  auroral  model,  see  Strickland  et  al}^). 

The  next  three  sections  present  our  full  sets  of  energy  loss  cross  sections  for  Nj, 
O2,  and  O,  including  references  to  all  individual  set  members.  As  noted  above,  total  cross 
sections  are  constructed  and  compared  to  electron  transmission  data.  Loss  functions  are 
also  constructed  and  compared  to  the  Bethe  formula.  A  discussion  section  (Sec.  5) 
completes  the  paper. 
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2.  ENERGY  LOSS  CROSS  SECTIONS  FOR  Nj 

Table  1  identifies  states  or  energy  loss  channels  corresponding  to  individual 
inelastic  cross  sections  of  N2.  For  each  entry,  the  table  also  includes  the  energy  threshold, 
location  of  the  cross  section  maximum,  the  value  of  this  maximum,  percent  contributing 
to  dissociation,  and  the  source  of  the  cross  section.  Tabulated  values  of  these  cross 
sections  (as  well  as  those  to  follow  in  Tables  2  and  3)  are  given  in  the  Appendices. 
Figure  1  shows  examples  of  measured  energy  loss  spectra  (from  two  separate 
measurements  as  noted  below)  with  an  energy  level  diagram  above  for  states  with  energy 
thresholds  within  the  illustrated  loss  region  (6  to  14.6  eV).  Such  data  are  obtained  by 
starting  with  a  beam  of  electrons  at  a  single  energy  and  measuring  its  energy  spectrum  at 
a  given  scattering  angle  after  passing  through  a  given  amount  of  N2.  High  lying  states 
from  Table  1  are  not  included  since  their  thresholds  (including  ionization)  lie  above  15 
eV.  Each  of  the  horizontal  line  segments  in  the  lower  portion  of  the  figure  is  identified 
with  a  given  state  and  shows  the  excitation  threshold  along  with  some  indication  of  the 
effective  range  of  energy  loss.  Unlike  the  localized  nature  of  energy  loss  for  an  atomic 
state,  here  loss  extends  to  several  eV  above  threshold  due  to  the  ability  of  an  impacting 
electron  to  leave  N^  in  one  of  several  vibrational  levels  of  a  given  electronic  state.  The 
loss  spectrum  below  12  eV  is  from  S.  Trajmar  (private  communication,  1995)  and  was 
obtained  for  electrons  with  an  incident  energy  of  40  eV  observed  at  a  scattering  angle  of 
20°.  The  vertical  scale  is  arbitrary  since  the  purpose  of  showing  the  data  is  to  simply 
illustrate  energy  loss  structure.  The  wings  of  the  Trajmar  loss  spectrum  have  been 
multiplied  by  10  to  show  the  structure  associated  with  the  various  triplet  states.  The 
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dominant  loss  for  the  given  incident  energy  is  seen  to  be  by  the  a'llg  state  responsible  for 
the  Lyman-Birge-Hopfield  band  system.  The  loss  spectrum  above  12  eV  is  from  Ratliff 
et  alP  for  100  eV  electrons  scattering  through  an  angle  of  15°.  The  structure  above  12.5 
eV  is  dominated  by  loss  to  numerous  vibrational  levels  of  the  b,  b\  c,  and  c’  states.  The 
first  figure  in  Ratliff  et  al?^  labels  the  peaks  by  vibrational  level.  The  magnitude  of  the 
Ratliff  spectrum  is  arbitrary  and  thus  no  significance  is  to  be  placed  on  its  strength 
relative  to  the  Trajmar  spectrum.  Cross  sections  are  obtained  from  data  such  as  these  by 
integrating  calibrated  spectra  over  angle  and  energy  loss  (see  papers  such  as  Ratliff  et 
al?^  and  Doering  and  Vaughan^^  for  more  information  on  the  derivation  of  cross  sections 
from  energy  loss  data). 

Figure  2  shows  a  comparison  between  the  total  cross  section  based  on  Table  1  and 
transmission  measurements  taken  from  the  review  by  Itikawa  et  af.  The  total  from  our 
work  is  comprised  of  the  three  curves  labeled  elastic,  total  ionization,  and  total  excitation. 
The  elastic  cross  section  was  also  taken  from  Itikawa  et  af.  The  total  and  elastic  cross 
sections  of  Itikawa  et  alf  are  based  on  available  measurements  with  adjustments  to 
account  for  offsets  among  the  various  data  sets.  References  to  original  data  may  be  seen 
in  the  paper  of  Itikawa  et  al.^  (See  also  Shyn  and  Carignan^®,  who  measured  the  total 
elastic  cross  section  from  1.5  to  400  eV.  This  reference  is  missing  in  Itikawa  et  alf). 
The  ionization  cross  section  comes  from  Rapp  and  Englander-Golden*'^  that  has  served  as 
the  standard  for  modeling  Nj  ionization  in  the  upper  atmosphere.  Recent  measurements 
by  Krishnakumar  and  Srivastava*^  are  in  close  agreement  with  Rapp  and  Englander- 
Golden*'*  (see  also  for  comparison,  the  derived  total  ionization  cross  section  by  Shyn^* 
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from  measurements  of  secondary  electrons  for  primary  electron  energies  from  50  to  400 
eV). 


The  excitation  cross  section  is  comprised  of  the  20  non-ionizing  components  in 
Table  1.  The  next  few  figures  present  these  components.  We  start  with  the  triplet  state 
cross  sections  in  Figure  3.  The  source  of  these  cross  sections  is  Cartwright  et  alP  with 
the  exception  of  that  for  the  C  state  which  is  Ajello  et  al?^  As  an  added  note,  Brunger 
and  Teubner^"*  have  recently  performed  energy  loss  measurements  similar  to  Cartwright  et 
al.  but  their  results  are  inconclusive  with  regard  to  integrated  cross  section  values  due  to 
the  restricted  range  of  scattering  angles.  Singlet  state  cross  sections  are  shown  in  Figure  4 
that  include  the  total  vibrational  cross  section  and  one  for  high-lying  states.  Figure  5 
shows  the  terms  comprising  this  latter  cross  section,  which,  from  Table  1,  are  seen  to 
come  from  Zipf  and  McLaughlin^^ 

In  Figure  6  we  show  a  comparison  between  the  total  dissociation  cross  section  that 
we  derive  from  our  full  set  of  cross  section  data  (see  Table  1)  and  those  measured  by 
Winters^^  and  most  recently  by  Cosby^®.  The  cross  section  data  of  Winters^’  have  been 
corrected  for  dissociative  ionization  by  using  the  recommended  values  of  dissociative 
ionization  cross  sections  from  Itikawa  et  af.  Although  the  values  of  both  measured  cross 
sections  for  energies  greater  than  15  eV  are  within  the  stated  error  limits,  the  systematic 
differences  seen  in  the  figure  may  be  due  to  an  additional  error  from  a  correction  for 
dissociative  ionization.  Our  cross  section  is  seen  to  be  in  overall  good  agreement  with 
both  measurements. 
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A  final  figure  before  discussing  Oj  addresses  energy  loss.  The  loss  function  based 
on  Table  1  is  compared  with  the  Bethe  formula  in  Figure  7  (see  Strickland  et  al}^\  Eq.  A9 
for  the  form  of  the  Bethe  expression).  The  comparison  is  restricted  to  energies  above  100 
eV  since  the  formula  begins  to  lose  its  validity  at  lower  energies.  The  cross  section  based 
loss  function  is  calculated  from: 


i(£)=X'no.(0+toi«»+l 
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Terms  are  as  follows: 
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The  three  components  in  Eq.  1  are  shown  in  Fig.  7.  The  importance  of  the  secondary 
electron  component  with  increasing  energy  is  a  reflection  of  the  increase  in  the  average 
energy  of  a  secondary  electron  as  the  energy  E  of  the  incident  electron  increases.  The 
differential  form  of  the  ionization  cross  section  is  given  by  Eq.  A4  in  Strickland  et  al}^ 

with  the  adjustable  parameter  E=  13  eV.  The  same  value  has  been  used  for  O2  and  O. 
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Excellent  agreement  is  achieved  with  the  Bethe  loss  function  which  gives  a  strong 
indication  that  ionization  is  being  conoctly  described  by  its  total  and  differential  forms  of 
the  cross  section.  While  excitation  dominates  the  loss  function  below  30  eV,  its 
contribution  at  higher  energies  falls  below  20%  where  comparisons  with  the  Bethe  loss 
function  become  valid.  Thus,  the  test  on  cross  sections  using  the  Bethe  formula  only 
weakly  addresses  excitation. 
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3.  ENERGY  LOSS  CROSS  SECTIONS  FOR  O2 


Table  2  shows  information  similar  to  Table  1  except  for  O2.  A  column  has  not 
been  included  for  contributions  to  dissociation  since  this  is  of  less  interest  than  for  N2 
(this  is  not  to  say  that  the  process  in  not  represented  within  the  full  set  of  cross  sections 

4 

presented  in  the  table).  Greater  interest  in  N2  is  due  to  chemistry  modeling  of  N(  S), 

2  2 

N(  D),  N(  P),  and  NO  within  the  AURIC  model  for  which  production  of  N  in  the  above 
states  by  dissociation  (by  photoelectrons  and  solar  photons)  must  be  specified.  Similar 
modeling  of  O  is  not  performed  due  to  O  being  one  of  the  dominant  species  in  the 
thermosphere  for  which  its  specification  is  given  by  a  model  such  as  MSIS-90  .  Like  Nj, 
all  excitation  cross  sections  listed  in  the  table  are  based  on  energy  loss  measurements 
except  for  those  belonging  to  Rydberg  states.  Unfortunately,  there  are  no  measurements 
available  for  these  states  and  consequently  we  have  adopted  a  theoretical  cross  section 
representing  total  Rydberg  excitation  as  compiled  by  Oran  and  Strickland^*^. 

Figure  8  shows  O2  energy  loss  spectra  similar  to  those  for  Nj  in  Fig.  1.  The 
spectrum  below  2  eV  is  from  Shyn  and  Sweeney"*®  for  10  eV  electrons  scattered  through 
an  angle  of  96°.  The  spectrum  above  7  eV  is  from  Shyn  et  a/."**  for  20  eV  electrons 
scattered  through  an  angle  of  156°.  Similar  to  the  spectra  for  N2,  the  results  in  Fig.  8  have 
arbitrary  scales  and  thus  there  is  no  significance  in  the  strength  of  the  low  energy  portion 
compared  to  that  above  7  eV. 
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Figure  9  shows  the  measured  total  cross  section,  the  corresponding  cross  section 
based  on  Table  2,  and  its  three  components.  The  measured  total  and  elastic  scattering 
cross  sections  come  from  the  most  recent  review  of  Kanik  et  al}^  (see  also  Sullivan  et 
Shyn  and  Shaip'^^  and  Wakiya'’^  for  original  studies  and  other  references).  Similar 
to  N2,  they  are  based  on  several  sets  of  measurements  with  adjustments  for  data  offsets 
among  these  sets.  The  total  ionization  cross  section  also  comes  from  Kanik  et  al}^  and 
includes  dissociative  ionization.  These  authors  compared  available  measurements  and 
recommend  those  recently  made  by  Krishnakumar  and  Srivastava*^  which  are  similar  to 
the  Rapp  and  Englander-Golden*'^  cross  section  within  the  stated  uncertainty  limits. 
Similar  results  are  also  derived  by  Shyn  and  Sharp'*^  from  secondary  electron 
measurements  from  threshold  to  300  eV.  The  excitation  cross  section  is  the  sum  of  the 
thirteen  from  Table  2.  Figures  10  and  11  show  the  individual  excitation  cross  sections 
where  a  sum  has  been  performed  over  the  vibrational  cross  sections.  Schumann-Runge 
dissociation  is  a  well  known  process  for  O2.  The  responsible  states  are  1  11^  and  B  2^ 
along  with  the  8.9  eV  channel  listed  in  Table  2 

Figure  12  presents  a  comparison  between  the  loss  function  calculated  with  eq.l 
and  the  Bethe  loss  function.  The  cross-section-based  loss  function  is  about  10%  below 
the  Bethe  loss  function  and  argues  for  an  increase  by  this  amount  in  the  magnitude  of  the 
ionization  cross  section.  As  with  N2,  energy  loss  by  excitation  is  minor  where  the 
comparison  is  being  made.  In  fact,  excitation  plays  a  weaker  role  in  O2  compared  to  N2. 
As  a  concluding  comment  in  this  section,  larger  errors  in  O2  energy  loss  cross  sections 


can  be  tolerated  in  photoelectron  and  auroral  electron  energy  loss  calculations  compared 
to  Nj  given  the  fact  that  there  is  much  less  O2  in  the  thermosphere. 
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4.  ENERGY  LOSS  CROSS  SECTIONS  FOR  O 


Table  3  shows  information  for  O  similar  to  that  in  Tables  1  and  2  for  N2  and  O2. 
There  are  many  high  lying  states  of  O  in  addition  to  those  explicitly  listed  (see,  e.g..  Fig. 
1  in  Laher  and  Gilmore^).  None  of  them  individually  accounts  for  significant  energy  loss 
based  on  the  measurements  of  the  investigators  referenced  in  the  table.  The  Rydberg 
cross  sections  represent  these  many  states  not  explicitly  accounted  for.  Figure  13  shows 
energy  loss  spectra  for  O  similar  to  previous  results  in  Figs.  1  and  8  for  N2  and  O2, 
respectively.  An  important  difference,  however,  is  the  local  nature  of  energy  loss  for  a 
given  electronic  state  compared  to  states  of  N2  and  O2  that  can  extend  over  several  eV 
due  to  vibrational  excitation.  The  spectrum  below  5  eV  is  from  Doering  and  Gulcicek"*^ 
for  30  eV  electrons  scattered  through  an  angle  of  120°.  The  spectrum  at  higher  energies 
is  from  Doering  and  Vanghn^^  for  100  eV  electrons  scattered  through  an  angle  of  4°. 
Again,  as  was  the  case  for  Nj  and  O2,  the  scales  for  the  two  spectra  are  arbitrary.  We 
have  included  energy  loss  data  in  Fig.  13  above  the  ionization  threshold  (13.6  eV).  The 
slow  rise  in  the  underlying  continuum  above  the  threshold  is  due  to  ionization.  The 
feature  at  9. 15  eV  has  been  added  based  on  data  from  Doering  and  Gulcicek^^  for  13.9  eV 
electrons  scattered  through  an  angle  of  50°.  Otherwise,  this  loss  feature  would  not  be 
discernible  given  the  incident  energy  (100  eV)  associated  with  rest  of  the  spectrum  above 
9eV. 


The  total  cross  section  for  O  and  its  components  are  shown  in  Figure  14.  The 
components  are  from  Itikawa  and  Ichimura*  for  elastic  scattering.  Brook  et  al}^  (also  see 
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Itikawa  and  Ichimura®)  for  ionization,  and  from  the  references  in  Table  3  for  excitation. 
Similar  to  N2  and  O2,  elastic  scattering  and  ionization  dominate  the  total  and  thus  the 
comparison  in  the  figure  demonstrates  little  about  the  accuracy  of  the  total  excitation 
cross  section.  The  next  three  Figures  (15-17)  show  the  excitation  cross  sections  listed  in 
Table  3.  Extrapolations  to  higher  energies  beyond  those  measured  are  based  on  expected 
fall-off  for  the  given  transitions. 

Similar  for  O2,  Rydberg  cross  sections  for  O  are  not  well  characterized 
experimentally,  and  consequently,  we  again  rely  on  theoretical  values  as  compiled  by 
Oran  and  Strickland^®.  The  Rydberg  cross  section  in  Figure  15  comes  from  Oran  and 
Strickland^®.  We  are,  however,  able  to  compare  with  other  theoretical  values  produced  by 
Laher  and  Gilmore^  who  have  carried  out  a  critical  review  of  inelastic  cross  sections  for 
O.  They  considered  more  than  sixty  individual  cross  sections,  including  nine  allowed  and 

+  4  +2 

twenty  nine  forbidden  Rydberg  series  cross  sections  with  transitions  to  O  (  Sq),  O  (  Dq), 
+  2 

and  O  (  Pq)  ion  cores  of  the  excited  states.  Since  no  cross  section  measurements  were 
available,  they  used  a  semi-empirical  formula  based  on  the  work  of  Jackman  et  alf'^  to 
estimate  individual  Rydberg  cross  sections.  Their  sum  is  included  in  Fig.  15  for 
comparison  with  our  representation.  Good  agreement  exists  below  30  eV  while  the  Laher 
and  Gilmore^  values  are  as  much  as  a  factor  of  two  higher  between  40  and  150  eV.  There 
are  large  uncertainties  associated  with  either  cross  section  and  we  thus  do  not  speculate  as 
to  which  is  more  accurate  based  on  the  respective  modeling  techniques.  Nevertheless, 
using  the  Laher  and  Gilmore^  cross  section,  our  total  cross  section  then  exceeds  the 
measured  one  by  about  10%  above  30  eV  whereas,  agreement  to  a  few  percent  is 
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achieved  using  our  representation.  This  provides  an  argument,  although  perhaps  weak, 
for  implementing  a  smaller  Rydberg  cross  section  above  30  eV  than  derived  from  the 
many  individual  terms  in  Laher  and  Gilmore^. 

Our  calculated  loss  function  for  O  is  shown  in  Figure  18  along  with  that  based  on 
the  Bethe  formula.  Reasonably  good  agreement  is  achieved  with  the  Bethe  formula 
although  the  differences  in  shape  and  magnitude  above  several  hundred  eV  suggest 
increasing  the  total  ionization  cross  section  by  perhaps  20%  and  decreasing  the  average 
energy  of  the  secondary  electron  per  collision  by  a  similar  amount. 
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5.  DISCUSSION 


From  the  information  presented  in  earlier  sections,  it  seems  clear  that  there  is 
much  less  uncertainty  in  the  ionization,  elastic,  and  total  cross  sections  of  a  particular 
species  compared  to  excitation.  It  is  difficult  to  assign  error  bars  due  to  the  many  sources 
of  the  information  although  our  own  assessment  is  that  15%  or  less  uncertainty  can 
probably  be  assigned  to  the  former  sets  and  more  than  30%  should  be  assigned  to  the 
latter  (excitation).  The  largest  source  of  uncertainty  for  excitation  is  in  the  Rydberg  cross 
sections.  The  situation  appears  to  be  most  satisfactory  for  N2  among  the  three  species 
being  addressed.  As  noted  earlier,  we  have  used  calculated  values  of  Rydberg  cross 
sections  for  O2  and  O  given  the  paucity  of  measured  values.  One  approach  to  specifying 
the  total  excitation  cross  section  is  to  1)  subtract  the  ionization  and  elastic  scattering  cross 
sections  from  the  total  cross  section  or  2)  subtract  the  ionization  cross  section  from  the 
total  inelastic  cross  section  (if  available).  We  do  not  recommend  such  subtractions  since 
the  calculations  involve  the  differences  of  similar  quantities  which  demand  greater 
accuracy  in  these  quantities  than  can  be  expected  at  this  time.  An  alternative  approach 
and  the  one  taken  in  this  work  is  to  compile  available  energy  loss  measurements  for 
specific  states  or  loss  channels,  supplement  them  where  necessary  with  theoretical  cross 
sections,  sum  this  total  set  along  with  the  ionization  and  elastic  scattering  cross  sections, 
and  compare  with  the  total  obtained  from  transmission  measurements.  While  such  an 
approach  does  not  assure  an  accurate  description  for  excitation  assuming  good  agreement 
between  totals,  it  is  nevertheless  a  worthwhile  exercise  for  assigning  some  degree  of 
confidence  to  the  overall  magnitude  of  the  total  excitation  cross  section. 
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An  important  aspect  of  our  N2  work  presented  in  Sec.  2  was  specifying  a  total 
dissociation  cross  section.  Based  on  its  magnitude  (see  Figs.  5  and  6),  dissociation 
accounts  for  approximately  80%  of  the  energy  loss  in  excitation  channels  for  electron 
energies  above  30  eV.  Given  the  strength  of  this  loss  channel  and  the  importance  of  odd 
nitrogen  to  a  number  of  aeronomical  problems,  we  provide  further  details  here  to 
supplement  the  discussion  in  Sec.  2.  Zipf  and  Mclaughlin^^  identified  nine  energy  loss 
channels  that  provide  most  of  the  contribution  to  the  N2  dissociation  cross  section  (see 

Table  1).  The  most  important  of  these  are  the  high-lying  states  and  the  family  of 
states  (terms  used  by  Zipf  and  Mclaughlin).  Cross  section  values  at  100  eV  have  recently 
been  measured  for  two  members  of  this  family.  Ratliff  et  al?^  addressed  the  state 
and  obtained  a  cross  section  equal  in  magnitude  to  29%  of  the  value  for  the  family  at  100 
eV.  James  et  al.*^  reported  a  value  for  the  state  equal  to  31%.  Since  these  new 

measurements  are  available,  we  have  constructed  cross  sections  for  these  states  using  the 
shape  of  the  cross  section  for  the  family  along  with  the  reported  magnitudes  at  100  eV. 
Having  removed  these  components  from  the  Zipf  and  McLaughlin^®  cross  section,  a 
residual  cross  section  for  the  remaining  members  was  constructed  simply  by  scaling  down 
the  Zipf  and  McLaughlin^®  cross  section  by  0.40.  While  the  decomposition  into  three 
new  members  with  the  same  shape  does  not  affect  photoelectron  or  auroral  electron 
energy  loss  calculations,  it  was  done,  nevertheless,  in  anticipation  of  further 
measurements  of  one  or  more  of  the  states  that  may  lead  to  some  differences  in 
shapes  among  the  components. 
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In  the  Introduction,  we  noted  that  several  investigators  have  compiled  sets  of 
energy  loss  cross  section  for  the  purpose  of  calculating  photoelectron  and  auroral  electron 
fluxes  and  that  differences  exist  among  these  sets.  Generally,  there  is  little  difference  in 
ionization  and  elastic  scattering  cross  sections  among  the  sets.  The  differences  occur 
among  the  excitation  cross  sections.  Examples  of  sets  that  have  noteworthy  differences 
are  those  of  Richards  and  Torr^'*’  and  Strickland  and  colleagues  (e.g.,  Strickland  and 
Meier^ ;  Strickland  and  Anderson^®).  Specifically,  the  total  N2  excitation  cross  section  of 
Richards  and  Torr^^  is  approximately  a  factor  of  two  smaller  above  15  eV  in  comparison 
to  that  of  Strickland  and  colleagues.  Richards  and  Torr^^  use  a  total  Nj  excitation  cross 
section  obtained  by  subtracting  a  total  ionization  cross  section  from  the  total  inelastic 
cross  section  of  Phelps  (private  communication)  as  appears  in  Stamnes  and  Rees^^  (more 
will  be  said  about  this  in  the  next  paragraph).  Strickland  and  colleagues  have  relied 
heavily  on  measurements  of  cross  sections  for  specific  states  or  specific  energy  loss 
channels  (e.g.,  Zipf  and  McLaughlin^^  ;  Cartwright  et  alP  ;  Ajello  and  Shemansky^'). 
Calculated  photoelectron  fluxes  above  15  eV  and  below  the  region  where  ionization 
begins  to  dominate  (above  about  50  eV)  by  Richards  and  Torr^'^’^^  are  approximately 
twice  those  of  Strickland  and  colleagues  due  to  these  differences.  This  has  led  to 
discussions  in  numerous  papers  about  the  accuracy  of  the  satellite  measured 
photoelectron  fluxes  by  Doering  and  colleagues  (e.g.,  Lee  et  Richards  and  Torr 

argue  that  the  fluxes  have  the  correct  magnitude  while  Strickland  and  colleagues  argue 
that  the  fluxes  should  be  reduced  by  a  factor  between  1.5  and  2.0  (e.g.,  Strickland  and 
Anderson^®).  Conway^^  presents  independent  information  that  also  argues  for  a  reduction 
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in  the  measured  photoelectron  fluxes.  A  resolution  to  the  problem  has  not  been  obtained 
as  of  this  writing.  It  is  important  to  note  that  similar  resuKs  (within  10%)  are  obtained  by 
us  using  either  the  cross  sections  in  this  paper  or  our  earlier  sets. 

Reference  was  made  in  the  above  paragraph  to  the  total  excitation  cross  section  of 
Phelps  in  the  Stamnes  and  Rees^^  paper.  Phelps  (in  a  private  communication  to  Stamnes 
and  Rees)  provided  several  cross  sections  for  triplet  states  and  one  for  total  excitation  to 
singlet  states.  Most  of  the  contribution  to  the  total  comprising  these  cross  sections  comes 
from  an  analysis  of  swarm  data.  Since  investigators  such  as  Stamnes  and  Rees  and 
Richards  and  Torr^'*  have  used  cross  sections  of  Phelps  and  colleagues  based  on  swarm 
data,  a  brief  discussion  of  recent  papers  addressing  such  data  is  presented  here.  The  Nj 
excitation  cross  sections  appearing  in  Stamnes  and  Rees  are  the  same  as  those  published 
by  Pitchford  and  Phelps^"^  with  the  exception  of  the  total  singlet  cross  section.  The 
version  appearing  in  Stamnes  and  Rees  is  about  twice  as  large  as  in  Pitchford  and  Phelps. 
The  increase  is  explained  by  Phelps  and  Pitchford^^  who  added  selected  high  threshold 
singlet  cross  sections  from  Zipf  and  McLaughlin^^.  Richards  and  Torr^^,  as  noted  above, 
subtracted  a  total  ionization  cross  section  (from  Kieffer  and  Dunn^^)  from  the  total 
inelastic  cross  section  in  Stamnes  and  Rees  (based  on  the  ionization  and  excitation  cross 
sections  provided  by  Phelps)  to  obtain  a  total  excitation  cross  section.  The  result  is  a 
smaller  cross  section  above  the  ionization  threshold  than  would  be  obtained  by  adding  the 
excitation  components  of  Phelps  since  the  ionization  cross  section  used  by  Richards  and 
Torr  in  the  subtraction  is  larger  than  assumed  by  Phelps  (Rapp  and  Englander- 
Golden*'^).  In  addition  to  a  larger  singlet  cross  section  by  Phelps  and  Pitchford^^ 
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compared  to  their  1982  value,  the  latter  paper  also  gives  a  larger  total  triplet  cross  section 
due  to  increasing  the  C  state  cross  section  by  a  factor  of  two.  The  more  recent  work  of 
Jelenkovic  and  Phelps^^,  based  on  stronger  electric  field  swarm  data,  retains  the  total 
singlet  cross  section  from  the  1985  work,  but  returns  to  a  total  triplet  cross  similar  to  that 
in  the  1982  work.  Compared  to  the  compilation  in  this  work,  the  total  singlet  cross 
section  is  -25%  higher  while  the  total  triplet  cross  section  is  -35%  smaller. 

Although  more  work  remains  to  be  done  on  characterizing  excitation  cross 
sections  for  N^,  O^,  and  O,  the  compilation  from  this  study  gives  a  comprehensive  set  of 
the  most  updated  energy  loss  cross  sections  that  should  serve  well  for  performing 
photoelectron  and  auroral  electron  energy  loss  calculations.  In  closing,  we  note  areas 
that  most  urgently  need  more  attention: 

1)  Measurements  of  the  total  scattering  cross  section  of  O  to  verify  and  understand  the 
limited  measurements  available  at  this  time.  Also  needed  is  an  experimental 
determination  of  the  O  elastic  scattering  cross  section. 

2)  Measurements  to  better  quantify  the  Rydberg  states  of  and  O. 

3)  Additional  selected  measurements  of  the  excitation  cross  sections  of  the  first  ten  states 

of  N,  (A"Z:,  C'n„,  triplets  and  a%,  wX 

a"^Zg  singlets)  to  resolve  differences  between  recent^'*  and  eeirlier  measurements^^’^^’ 
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4)  Measurements  of  cross  sections  for  various  'n„  states  of  Nj  to  verify  the  measurement 
of  Zipf  and  McLaughlin^^  for  the  sum  of  cross  sections  for  these  states.  As  noted 
earlier,  single  energy  measurements  have  recently  been  made  for  the  b^Uu  and  c‘n„ 

states.  More  measurements  are  needed,  especially  near  30  eV  where  these  cross 
sections  peak. 
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TABLES 


Table  1.  Ionization  and  excitation  cross  sections  for  N^.  We  are  using  these  to  account 
for  all  energy  loss  to  electrons  impacting  on  in  calculations  of  photoelectron 
and  auroral  electron  fluxes. 

Table  2.  Ionization  and  excitation  cross  sections  for  Oj.  Similar  to  Table  1,  these  are 
intended  to  account  for  all  energy  loss  to  electrons  impacting  on  O^. 

Table  3.  Ionization  and  excitation  cross  sections  for  O.  Similar  to  Table  1,  these  are 
intended  to  account  for  all  energy  loss  to  electrons  impacting  on  O. 


FIGURES 

Figure  1.  Measured  energy  loss  spectra  for  and  energy  level  diagram  showing 
electronic  states  contributing  to  loss  in  the  given  energy  range.  See  text  for 
details  on  incident  electron  energies,  scattering  angles,  and  references. 

Figure  2.  Comparison  between  our  derived  and  measured  total  N2  cross  section.  Also 
shown  are  the  components  of  our  constructed  cross  section  that  include  elastic, 
total  ionization,  and  the  sum  of  excitation  cross  sections  from  Table  1. 

Figure  3.  Excitation  cross  sections  for  the  six  triplet  electronic  states  of  N2. 

Figure  4.  Excitation  cross  sections  for  the  ten  singlet  electronic  states  of  N2.  Vibrational 
cross  sections  are  summed  into  a  single  cross  section  whose  magnitude  is 
reduced  by  a  factor  of  10  to  be  fully  displayed  in  the  given  panel.  Also  shown 
in  the  lower  panel  is  the  sum  of  the  cross  sections  for  the  high  lying  states. 

Figure  5.  Individual  excitation  cross  sections  for  the  high  lying  states. 

Figure  6.  Derived  and  measured  cross  sections  for  dissociation  (excluding  dissociative 
ionization).  Error  bars  refer  to  the  Cosby  data. 

Figure  7.  Loss  function  from  this  work  compared  with  the  Bethe  formula.  Also  shown 
are  its  components  from  this  work. 

Figure  8.  Examples  of  energy  loss  spectra  for  and  identification  of  states  producing 
the  exhibited  structure.  See  text  for  details  on  incident  electron  energies, 
scattering  angles,  and  references. 

Figure  9.  Total  cross  section  information  similar  to  Figure  2  except  for  O^. 
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Figure  10.  Cross  sections  for  vibration  excitation  (summed  over  components  identified  in 
Table  2),  for  the  singlet  states,  and  for  a  combined  measurement  over  the 
i-A'^A^+c'Sy  states. 

Figure  11.  The  remaining  excitation  cross  sections  from  Table  2. 

Figure  12.  Loss  function  information  similar  to  that  in  Figure  7  except  for  O^. 

Figure  13.  Examples  of  energy  loss  spectra  for  O  and  identification  of  states  producing 
the  exhibited  structure.  See  text  for  details  on  incident  electron  energies, 
scattering  angles,  and  references. 

Figure  14.  Total  cross  section  information  similar  to  Figure  2  except  for  O. 

Figure  15.  Cross  sections  for  singlet,  quintet  and  Rydberg  states  of  O.  Rydberg  cross 
sections  (as  summed  by  us)  from  Laher  and  Gimore^  are  also  shown  for 
comparison. 

Figure  16.  Cross  sections  for  the  triplet  states  of  O. 

Figure  17.  Autoionization  cross  sections  of  O  for  the  following  transitions: 

2s  V  3s''  and  4d'  ^P°. 

Figure  18.  Loss  function  information  similar  to  that  in  Figure  7  except  for  O. 
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Cross  section  (cm^  Cross  section  (cm*) 


FIGURE  3 
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FIGURE 


(jUio)  uopoas  SS0J3 
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Energy  (cV)  Energy  (eV) 

FIGURE  9  FIGURE  10 


(as)  X2j3U3 
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Rydbergs  (Ref.®) 
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Loss  Function  (eV  cm^) 


Energy  (eV) 

FIGURE  18 
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Table  1 


Excitation/ 

Ionization 

Threshold 

(eV) 

E  (eV) 

a  (10'*) 

cm^ 

%  Contribution 
to  Dissociation 

Reference 

N;  {Total) 

15.6 

100 

252 

— 

12, 14 

li^{X^*)vib  {Total) 

0.9 

2.0 

1540 

— 

6 

N,(a'np 

9.1 

18.0 

26.9 

12 

51, 59, 60 

14.2 

60.0 

13.6 

83 

61 

N,(c:'x) 

12.9 

80.0 

12.4 

15 

61 

N,  (!>‘n.) 

12.6 

40.0 

21.2 

96 

28, 36 

N,(c‘n) 

12.9 

40.0 

22.0 

100 

49, 36 

8.4 

15.0 

10.4 

— 

35 

N2(a"'z;) 

12.3 

20.0 

5.8 

— 

35 

N2(iv‘A„) 

8.89 

13.0 

11.7 

— 

35 

N,(A^:) 

6.2 

17.0 

22.0 

— 

35 

N,(5^n,) 

7.4 

11.5 

29.5 

— 

35 

N,(c"np 

11.0 

14.0 

36.0 

50 

33, 36,61 

N,(W*A^ 

7.5 

16.5 

38.0 

— 

35 

N,(5^:) 

8.0 

15.0 

12.5 

— 

35 

12.0 

24.0 

0.8 

— 

35, 36 

N,(i5.5cVpeaA:) 

16.4 

40.0 

25.0 

100 

36 

N,(VIA0 

23.7 

100 

15.5 

100 

36 

li^{17.3eVpeak) 

17.4 

40.0 

10.5 

100 

36 

Nj  {N  Ryd  atoms) 

40.0 

88.0 

3.4 

100 

36 

Nj  (A/j  triplet  manifold) 

11.0 

17.0 

12.0 

100 

36 

Other  'riu  states 

12.6 

40.0 

30.0 

— 

36 
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Table  2 


Excitation/ 

Ionization 

Threshold  (eV) 

E  (eV) 

a_  (lO  ")  cm^ 

Reference 

O2  {Total) 

12.1 

100 

291 

10, 13, 14 

02(X’z;)v=l 

0.3 

8.8 

35.1 

40 

02(Z^pv  =  2 

0.4 

9.9 

17.6 

40 

0,(X^;)v=3 

0.6 

9.7 

7.7 

40 

02(X^;)v  =  4 

0.8 

9.3 

4.5 

40 

7.6 

20.0 

6.8 

41 

0,{^:)SR 

8.3 

19.3 

60.6 

41 

0^  (8.9  eV peak)  SR 

8.9 

24.3 

14.5 

41 

Oj  (second  band)  (^2) 

10.3 

21.8 

1.0 

65 

oM\) 

1.0 

7.4 

10.6 

63,64 

1.6 

6.9 

3.5 

63,64 

Oj  (longest  band)  (^2) 

10.0 

24.4 

6.8 

65 

02(a"2:+a'X+c*5:;) 

4.5 

8.0 

15.3 

7,64 

Oj  (Rydbergs) 

16.0 

32.0 

140.0 

18 

Table  3 

Excitation/ 

Ionization 

Threshold 

(eV) 

E  (eV) 

Reference 

0*  (Total) 

13.6 

100 

138 

8, 15 

2s^  2p^  '£) 

2.0 

6.0 

28.3 

47,69 

2s^2p^  ‘5 

4.2 

9.0 

3.3 

47,70 

2^2p^3s^S° 

9.3 

13.7 

8.9 

46 

2s'^2p^3s^S° 

9.5 

29.5 

9.1 

66 

2s^2p^3p^P 

10.7 

16.0 

3.3 

67 

2s^2p^3p^P 

11.0 

24.0 

12.4 

67 

2^2p^3d^D° 

12.1 

45.6 

36.6 

66 

2s^2p^3s'^D 

12.5 

49.0 

5.9 

68 

2s^2p^4d^D° 

12.8 

45.1 

2.0 

68 

2s'^2p^5d^D° 

13.0 

47.5 

1.0 

68 

2^2p^Ad^P° 

16.0 

45.1 

2.8 

68 

2s^2p^  3po 

15.0 

45.1 

13.2 

68 

2?  2/  3f" 

14.0 

44.8 

13.3 

68 

Rydbergs 

14.0 

25.0 

30.0 

18 
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Table  A1 


Ionization 

Dissociation 

Vibrational 

E 

a 

E 

a 

E 

c 

(eV) 

(10'”  cm^ 

(eV) 

(10'”  cm^ 

(eV) 

(10"”  cm^ 

— 

11 

.070 

— 

— 

12 

.328 

16 

.113 

14 

1.77 

18 

1.33 

16 

3.32 

.111 

20 

2.68 

20 

6.59 

8.84 

30 

9.78 

30 

12.08 

101. 

50 

18.82 

50 

12.84 

135. 

100 

25.20 

100 

11.43 

3.0 

58.6 

200 

22.58 

200 

8.48 

4.0 

13.7 

500 

14.55 

500 

5.03 

5.0 

2.34 

1000 

9.21 

1000 

3.16 

6.0 

.366 

Table  A2 
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Table  A3 


Table  A4 


c 

'n„ 

g 

Other  states 

E 

a 

E 

a 

E 

0 

E 

a 

(eV) 

(10'"  cm') 

(eV) 

(10'"  cm') 

(eV) 

O 

o 

B 

(eV) 

(10'"  cm') 

—— 

11 

.507 

13 

.031 

— 

— — 

13 

,002 

12 

.688 

14 

.182 

13 

.002 

14 

.163 

14 

.910 

16 

.368 

14 

.216 

16 

.549 

16 

.824 

20 

.551 

16 

.724 

20 

1.26 

20 

.558 

25 

.426 

20 

1.66 

30 

2.17 

30 

.334 

30 

.304 

30 

2.86 

50 

2.13 

50 

.193 

50 

.152 

50 

2.81 

100 

1.70 

100 

.096 

100 

.069 

100  . 

2.25 

200 

1.21 

200 

.045 

200 

.033 

200 

1.59 

500 

.705 

500 

.016 

500 

.013 

500 

.930 

1000 

.458 

1000 

.008 

1000 

.007 

1000 

.605 
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Table  B1 


Vibrational 


E 

(eV) 

a 

(10"'^  cm' 

1 

.145 

2 

.381 

3 

.720 

4 

1.17 

6 

3.24 

9 

5.98 

12 

2.59 

15 

.755 

20 

.191 

30 

.030 

Rydbergs 


17 

.281 

20 

4.18 

30 

14.09 

50 

12.36 

100 

7.59 

200 

4.34 

500 

2.18 

1000 

1.26 

Ionization 


E 

a 

(eV) 

(10*'’ cm' 

13 

.200 

16 

1.09 

20 

2.93 

30 

8.08 

50 

18.24 

100 

29.08 

200 

27.46 

500 

17.10 

1000 

10.54 

Table  B2 
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0. 

*n 

'!!■ 

. 

« 

s 

b  e" 
o 

y*N 

w  ^ 

18  .019 

20  .038 

25  .085 

30  .140 

50  .264 

100  .178 

200  .081 

500  .021 

800  .010 

ft, 

X 

% 

cn 

£ 

O  e" 
o 

m  > 

^2. 
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